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Abstract In this review we describe recent observational and theoretical de¬ 
velopments in our understanding of pulsar winds and pulsar-wind nebulae 
(PWNe). We put special emphasis on the results from observations of well- 
characterized PWNe of various types (e.g., torus-jet and bowshock-tail), the 
most recent MHD modeling efforts, and the status of the flaring Crab PWN 
puzzle. 
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1 Observations of PWNe. 

1.1 Introduction 

Only ^ 1% of the total pulsar spin-down luminosity is emitted as pulsed 
electromagnetic radiation, the majority of the spin-down luminosity of a pulsar 
being carried away by a relativistic and highly magnetized pair plasma. These 
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particles are generally believed to be accelerated and randomized in their pitch 
angle either upstream or at the pulsar wind termination shock. The radiation 
produced but these particles downstream of the termination shock is often 
seen as a pulsar-wind nebula (PWN). 

Most of the recent progress in our understanding of PWNe has been spurred 
by X-ray and TeV 7 -ray observations. The majority of PWNe has been dis¬ 


covered in one of these bands, and many are seen in both (see Kargaltsev et al 


2013b). To study PWN emission, it is important to disentangle the pulsar and 


pulsar wind contributions either by spatially resolving the nebula from the 
pulsar or by isolating the PWN component in the spectrum (e.g., the PWN 
contribution is expected to dominate in TeV). High-resolution images from 
Chandra X-ray Observatory {Chandra hereafter) revealed two dominant PWN 
morphologies: torus-jet and bowshock-taiQ In addition, a few objects with 
varying and puzzling morphologies can be seen in Chandra images ( [Kargalt^ 
sev and Pavlov||2QQ8 ). The PWN properties (size, morphology, and spectrum) 


can be expected to depend on the pulsar parameters (spin-down properties, 
pulsar velocity, and the angles between the spin, magnetic dipole, and velocity 
vectors) and on the environment (e.g., ambient pressure, magnetic field, and 
radiation field). The limited angular resolution of the ground-based TeV arrays 
(such as H.E.S.S. and VERITAS) does not allow us to detect TeV emission 
from the same particles that produce bright and compact X-ray nebulae in 
the vicinity of the pulsar. Instead, TeV images reveal much larger structures 
filled with the aged particles that may have accumulated over substantial part 
of pulsar’s lifetime (see de Jager and Djannati-Ataf| 2009). The TeV emis¬ 
sion is usually attributed to inverse Compton (IC) scattering of background 
optical/IR photons off aged electrons, although, in denser environments, the 
contribution of neutral pion decay to the 7 -rays emission could play an im¬ 
portant role (if the pulsar winds indeed contain the so far elusive relativistic 
protons). 

A general overview of PWN physics and X-ray observations was presented 
by Gaensler and Slane (2006) and Kargaltsev and Pavlov (2008), while PWN 
theory was recently reviewed by Amato (2014) and Bncciantini| ( [2014] ) . Here we 
focus on some of the most recent observational result^ and their implications 
(§1), discuss the latest theoretical advances in MHD modeling (§2), and review 
the non-MHD scenarios that can explain the puzzling Crab PWN flares (§3). 


1.2 The Crab and Vela PWNe as prototypes of young PWNe in SNRs. 

The Crab and Vela PWNe are often considered to be archetypal representatives 
of torus-jet PWNe. Since the Vela PWN is a factor of 10-20 older than the 


^ Such classification is possible only for sufficiently bright and relatively nearby PWNe. 

^ A more det ailed overview of t he obse rvational properties of population of relic PWNe 
can be found in [Kargaltsev et ~aL\ ( [20131^ . 
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Crab PWN, one can look for evolutionary changes by comparing the twc[^ 
PWNe of this type are usually found around young pulsars whose velocities 
are smaller than the speed of sound inside their host SNRs. These environments 


can be characterized by relatively high pressures and temperatures (Bamba 
et al.||2QlQ ). There is also evidence that pulsars powering torus-jet PWNe are 


likely to have substantial misalignment between the spin and magnetic dipole 
axes (e.g., Crab and Vela pulsars; Moffett and Hankins 1999 and Johnston 


et ah I 200^ respectively) which may play a pivotal role in formation of this 


type of morphology. 


1.2.1 Multiwavelength properties of the Crab 


The Crab nebula has been studied with nearly all major telescopes since its 
discovery (see Hester|[2QQ8 and references therein). However, it was not until 


the Hubble Space Telescope {HST; Hester et al.| 1995) and Chandra (Weis- 
skopf et al.||2QQQ| era when the intricate and complex structure of the nebula 


was revealed (see Figure for the feature nomenclature introduced by 


Hester 


These observations have also shown that the bright inner part of 


the nebula is very dynamic, with apparent velocities corresponding to up to 
0.5c (in projection onto the sky) as measured, e.g., from the shifts in wisp po¬ 
sitions ( [Hester et al.|2QQ'2 ). The changes in the nebula are more complex than 
simple translational motion (e.g. steady expansion). They include variations 
in brightness (e.g.. Inn er Knot; [Melatos et aT]|2QQ5 ) and shape (e.g.. Sprite; 
Bietenholz et al.|[2QQ4 Hester et al. 2QQ2| ). The wisp shapes can also be very 


different and while most of the wisps can be described as a ripple pattern with 
ripples moving away from the pulsar some of the wisps appear at the same 
location (e.g.. Thin Wisp in Figure [^. The prominent south-eastern (SE) jet 
(see Figure]^ also shows quite remarkable changes in its shape, based on X- 
ray images (taken over 14-year baseline), which could be explained by either 
precession of the curved jet or by the motion of kinks along the jet (Weisskopf 
2012 Deshpande and Radhakrishnan|2QQ7 ). Finally, the unexpected detection 


of 7 -ray flares by Fermi LAT and AGILE (Buehler et al. 2012 Abdo et al. 


2011b Tavani et al. 2011) suggests a significa nt energy release rate (which 


can reach 4 x 10^^ erg s Buehler et al. 2012) on timescales of hours (Mayer 


et al. 2013); however, it was not possible so far to pinpoint the location of 
the flaring region because of the lack of a “smoking gun” at lower frequencies. 
Consequently, the lack of information about the site of the flare in the PWN 
has led to a variety of models being suggested (see § . Given the lack of con¬ 

temporaneous variability at lower frequencies, it might be possible that some 
of the energy released during the process associated with the y-ray flares will 
manifest itself as more gradual flux changes at lower frequencie^ occurring on 


^ Since the PWN properties and evolution depend on the environment, one should not 
forget that the progenitor and SNR properties may be quite different for the Crab and Vela 
pulsars. 

See animation at http://home.gwu.edu/~kargaltsev/Crab.html 
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much longer timescales (e.g., hard X-ray variability reported by Kouzu et al 
2013 Wilson-Hodge et aL]|2Qll ). 


As a baseline for further comparison a multi wavelength (MW), high-resolution 
snap-shot of the PWN was obtained within a single day (Krassilchtchikov et 
al. 2015 in preparation; K-1-15 hereafter). Figure]^ shows Chandra^ HST (NIR, 
optical), and Karl G. Jansky Very Large Array (JVLA) images from this latest 
MW campaign. These data are ideally suited for measuring the contempora¬ 
neous spectra of the prominent PWN features. Both the radio and optical 
(broad-V band) images reveal prominent filaments of which only some are co¬ 
incident (which implies different MW spectra). Also, the SE jet (“1” in Fig. 
[^, which is prominent in the X-ray image, appears to have counterparts in 
the optical and NIR images but not in the radio image, in agreement with 
the earlier findings of [Bietenholz et ah (2004). On the other hand, a jet-like 
structure (“2” in Fig. located to the east of the SE X-ray jet in the JVLA 
image does not have an X-ray or NIR counterpart and coincides with one of 
the thermal optical filament^ Finally, in the NIR image, to the west of the 
bright X-ray jet (“1” in Fig.l^ there is another bright linear feature (“3” in 
Fig.§ extending southward from pulsar. The feature has a radio counterpart, 
a very faint X-ray counterpart, and does not coincide with any of the opti¬ 
cal filaments. The faintness of the X-ray counterpart suggest that emission 
from this feature is a synchrotron continuum may be produced by a cooled 
population of particles. It is, however, difficult to explain the presence of this 
low-energy feature in the conventional axisymmetric paradigm with equatorial 
and polar outflows, where only one jet and one counter-jet are expected^ On 
the other hand, the “line-free” (F550M) optical image shown in Figure [^sup¬ 
ports the axial outflow paradigm by revealing better than ever the other side 
of the axial “backbone” of the torus (dubbed as a counter jet, or NW jet). We 
note that at the brightness/contrast level chosen in the F550M image (Fig.[^ 
bottom left panel) the counterpart of the SE X-ray jet is barely seen, while 
the NW jet and the SE “jet-like” feature (see above) are clearly seen. Similar 
to the SE X-ray jet, the counter-jet appears to be hardly discernible in the 
radio images (see e.g., right panel in Fig. b ut it stands out in the 2-epoch 
difference image produced by Bietenholz et al. (2014) and shown in the bottom 
right panel of Figure Finally, the counter-jet is so faint in X-rays (if preset 
at all) that it cannot be discerned from the torus emission. Therefore, the 
frequency-dependent differences in brightness between the jet and counter-jet 
appear to be more complex than those expected for a simple scenario with the 
frequency-invariant Doppler boost (see e.g., Schweizer et al.|[2Q13 ). 


^ Overall, the large degree of correlation between the radio structure and optical filaments 
suggests than most of the radio emission is related to the SNR filaments. 

® However, one can imagine that the jet stayed at one position (“3” in Fig. [^ for a long 
time and then relatively quickly moved to the other position (“1” in Fig.j^ 
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K+15 found that the location of the bright optical/NIR wisp (“4” in Fig. 
only approximately coincides with the X-ray rin^ and the wisp brightness 
in the NIR image drops much faster with the distance from the symmetry 
axis of the nebula compared to the X-ray ring brightness. This supports the 


Schweizer et al. (2013) findings (based on earlier optical and X-ray monitoring) 


who concluded that the X-ray and optical emission must be produced by 


different populations of particles. Furthermore, according to Schweizer et al. 


(2013), the fits with the Doppler-boosted tilted ring model require noticeably 
different (higher) flow velocities for the optical wisps (~ 0.9c) compared to the 
X-ray wisps, which made Schweizer et al.|(|2013 ) question the simple “boosted- 
ring” model (see, however, Olmi et al.||2015 ). K+15 also found that for most 
individual features of the PWN (e.g., wisps) the NIR-optical-FUV spectra are 
harder than the contemporaneously measured X-ray spectra suggesting either 
an additional narrow spectral component or other kind of complex behavior 
between optical and soft X-rays. 

The overall morphology of the Crab nebula has been reproduced in the rel¬ 
ativistic magnetohydro dynamic simulations with anisotropic energy flux (most 
recently by Volpi et ^|2QQ8 ; Camus et al.||2QQ9 Forth et ar]|2Q14 ). The sim¬ 
ulated images display wisp-like features and other variable structures along 
the symmetry axis. In the MHD simulations these features are formed due 
to the combination of flow dynamics (i.e. fluid vortices produced near the TS 
due to magnetosonic oscillations; see e.g., (Olmi et al. 2015) and references 


therein) and Doppler boosting. The emission from the oblique termination 
shock that has been also associated with the Crab’s Inner Knot (Komissarov 
and Lyubarsky 2QQ3[ [Lyutikov et al. 2Q15[ Yuan and Blandfor^ 2Q15| see, 


however Melatos et al.|2QQ5 and references therein for alternative interpreta¬ 
tions). On the other hand, some challenges to the MHD models still remain 
(Biihler and Blandford 2014). The predicted bright arc due to the emission 


from the relativistic post-shock flow originating from the termination shock 
(Komissarov and Lyubarsky 2003 Del Zanna et al]|2006 ) is absent in the 
X-ray images where we instead see the patchy (likely consisting of multiple 
knots) inner ring which may or may not appear to be brighter on the north¬ 
western side (depending on the observation epoch). In this sense the optical 
(or NIR) images featuring a bright wisp NW of the pulsar appear to be in a 
better correspondence with the predictions of the MHD models. The models 
also predict co-spatial small-scale structures in the optical and X-rays while 
this is generally not observed (see above). Only parts of the X-ray inner ring 
are seen in the optical and most (if not all) optical wisps also do not appear 
to have co-spatial X-ray counterparts (e.g., the Thin Wisp[^ labeled in the 
top panel of Figure [T] is lacking any nearby counterpart while other wisps are 
offset from the possibly associated X-ray bright features; see [Schweizer et ahl 


^ The outer edge of the optical/NIR wisp is about further away from the pulsar, with 
the X-ray ring emission trailing behind or possibly being sandwiched between the bright 
wisp and the fainter wisp. 

® Here we are following nomenclature introduced by Hester et al. ( 1995| , see Fig. 0 (top 
panel). 
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Fig. 1 The top panel introduces the most prominent features of the Crab PWN and their 
“conventional” names (from ( [Hester et al.|199^ ). The bottom left panel shows the 12.5 ks 
exposure image obtained with HST ACS F550M. The image has been produced by a combin- 
ing series of a uxiliary images obtained during the 09/2005-12/2005 polarimetry campaign 
( |Hester||2008| >. The F550M filter avoids any strong emission lines and provides a relatively 
unobstructed view of the synchrotron nebula. The “?” mark enigmatic feature discussed in 
§1.2.1 (also labeled as “3” in Fig. top right panel). Notice that the feature labeled as 
the counter-jet in the F550M image also appears to show large variability in the difference 
image (shown in the bottom right panel) produced from 2 JVL A observations obtained on 
2001 April 16 and 2012 August 26 (from [Bietenholz et al.|2014| >. 
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Fig. 2 False color MW images of the Crab PWN (see the legends in the panels) based 
on the observations described and analyzed in K+15. Numbers refer to the PWN features 
mention in the text. 


2013 and K+15). Therefore, it yet remains to be shown whether more ad¬ 


vanced models can fully capture the rich MW structure and variability of the 
Crab PWN. It seems that a complex injection spectrum may be required to 
achieve this, hinting that there may be multiple acceleration sites throughout 
the PWN with possibly different acceleration mechanisms. The most recent, 
advanced 3D models predict somewhat disordered structure of the magnetic 
field and suggest the need for the in-situ particle acceleration outside the ter¬ 
mination shock region (Porth et al. 2014). Diffusion transport may become 
more important in the case of disordered magnetic field. 


1.2.2 Multiwavelength properties of the Vela PWN 

The Vela pulsar is a factor of 20 older than Crab, hence some evolutionary 
differences are expected. In addition, the differences can be attributed to a 
different progenitor type, different ISM, and different properties of the pulsar 
(e.g., magnetic field, angles between the spin and magnetic dipole axis, or 
the orientation and magnitude of the pulsar velocity). Figure 4 (bottom left) 
shows a deep Chandra ACIS image of the Vela PWN (see also zoomed-in view 
of the compact nebula in Fig. produced by combining images from the 
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Fig. 3 Left: 12.5 ks exposure image of the Crab PWN obtained with HST ACS F550M. 
The image has been produced by combining a series of auxiliary images obtained during the 
09/2005-12/2005 polarimetry campaign |Hester| ( |2008| > . The F550M filter avoids any strong 
emission lines and provides a rel atively unobstructed vi ew of the synchrotron nebula. Right: 
JVLA image of the Crab (from [Bietenholz et al.|2014| >. 


latest observational campaign (comprised of eight 40 ks observations taken 
with one-week intervals; see [Durant et al.|[2QT3l The overall morphology of 
the bright compact PWN can be described as an axisymmetric double-arc 
structure with two axial jets having different brightnesses and widths. The 
bright, compact X-ray PWN is located inside the larger double-lobe radio 
PWN (Dodson et al. 2QQ3| ). Interestingly, the radio lobes appear to be filled 
by fainter X-ray emission which is particularly clearly seen in the ACIS hard 
band (1-8 keV; see Figure 4) thus suggesting that the radiating particles have 


not cooled much. Therefore, the pure advection model proposed by Kennel 


and Coronitij ( |1984a| ) may need to be augmented with some other transport 
mechanism (e.g., diffusion) capable of moving the energetic particles away 
from the pulsar more rapidly. We also note that the bright, double-lobed radio 
nebula is surrounded by a much larger (^ 2° in diameter) radio-emitting 
structure called Vela X (see §1.5.4) which is filled with bright filaments and 
fainter diffuse continuum (Figure]^ bottom right). 


1.2.3 The Crab and Vela PWNe: similarities and differenees 

The double-arc X-ray morphology of the Vela PWN is quite different from the 
Crab PWN, with its single ring, and torus (likely comprised of multiple wisps) 
seen in the ACIS images (see Fig.|^. These differences can hardly be attributed 
to the larger age of the Vela PWN (because the X-ray emission from compact 






















Pulsar-Wind Nebulae 


9 



8‘'42'" 39*" ae*" sa*" 30"’ 27"" 

RA (J2000) 


Fig. 4 Comparison of the Crab and Vela PWNe. The top panels show the combined X- 
ray (red) and radio (blue) images of the Crab (top left) and Vela (top right) PWNe. The 
bottom panels show X-ray and radio images of the Vela PWN. The Chandra ACIS image 
(bottom left) shows that faint X-ray emission seems to fill in the radio lobes (shown by the 
conto urs). The larger ra dio image on the right shows the entire Vela X complex (radio image 
from ( [Frail et al.||1997| ) within the Vela SNR. The inset shows the compact radio nebula 
(blue color) and brightest part of the X-ray nebula (red). 


PWN mostly comes from freshly injected electrons) or different ambient (SNR) 
medium properties (these could become progressively more important further 
away from the termination shock). It is also unlikely to be due to the difference 
in the angles between the pulsar spin axis and the line of sight because these 
angles are believed to be similar (127° for Vela [Helfand et al.|2QQi] and 120° 
for Crab Weisskopf et al.||2Q12| ). The only other important parameters could 
be the angle between the spin and magnetic dipol^ axis (still rather poorly 


^ In principle, the magnetic field may deviate from the dipolar configuration more for the 
Vela pulsar than for the Crab pulsar. Indeed, braking indices, n, of the Crab (n = 2.5) and 
Vela (n = 1.4) are very different, and the Vela pulsar is much more “glitchy” compared to 
Crab. 
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constrained, e.g., ^ 43° for Vela Johnston et al. 2QQ5| and ~ 45° — 70° for Crab 
Lyne et al.|2Q13 ) and different pulsar velocitie^^^|(the projected onto the sky 
velocities are OOdo.a and ^ 120(^2 km s“^ for Vela and Crab, respect ivel}p^ . 
These velocities are smaller than the typical sound speed inside the young 
SNR, and therefore the Mach number A4 should be ^ 1 for both pulsars^ 
Thus the differences in the compact PWN morphologies are more likely to^ 
attributed to the different angles between the spin and magnetic dipole axis. It 
is also possible that the different degrees of deviation of the NS magnetic field 
from that of an ideal, centered dipole have some impact on the PWN. Even for 
the compact parts of the PWNs some of the differences in X-ray morphologies 
could still be attributed to the longer synchrotron cooling time for the Vela 
PWN which is expected to have weaker magnetic field ( Pavlov et al.||2QQ3 ). 

We also note that the physical connection between the bright inner jets in 
the Vela PWN and the fainter large-scale outer jets (see e.g.. Fig. [^bottom left 
panel) has not been established yet. The bright SE inner jet of the Vela pulsar 
suddenly appears out of the orthogonal bar-like feature (shock in the polar 
backflow? [Komissarov and Lyubarsky 2003) at about 5.4" from the pulsar 
and then nearly as abruptly fades away at about 10.5" from the pulsar (see 
Fig. left panel). In the very deep Vela PWN image the outer jets are visible 
up to much larger scales 2 — 3'); however, we do not see any smooth 
transition from the inner axial jets and hence we cannot establish a firm link 
between the two (except that both structures are extending along the PWN 


The parameter that determines to what degree a PWN is affected by the pulsar motion 
is the ratio of the pulsar velocity to the local ambient sound speed (Mach number At). The 
medium within the younger Crab SNR is hotter than in the Vela SNR. 

Here the distances are scaled as do .3 = d/(300 pc) and d 2 = d/2,000 pc. 

It is likely that the M ach number is som ewhat larger for the Vela PWN where we see 
some effect of the motion ( [Pavlov et al.|2Q03| >. 
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sym metry axij^. Not e that simulated X-ray images based on MHD models 
(e.g. Del Zanna et al.|[2QQ6 [Camus et al.||2QQ9 ) do not show such structures. 
In the Crab PWN the dynamical feature called “Sprite” (possibly an analog 
of the bar feature in the Vela PWN) seems to be the place from which the 
SE jet originates. However, the Crab’s SE jet does not undergo any dramatic 
transitions in brightness until it bends and terminates at about 1.2' from the 
pulsar. Morphologically, this “kinked” jet resembles the outer jets of the Vela 
PWN rather than the straight and bright inner jets. This leaves no obvious 
analog to the bight inner jets of the Vela PWN in the Crab PWN (see, however, 
Deshpande and Radhakrishnan|2QQ7 ). As it has been discussed by [Pavlov et al.| 
"(2003) and Durant et al. (2013), the fact that the NW outer jet of the Vela 


PWN is brighter than the SE jet is at odds with the 3D orientation inferred 
from the arc brightness distribution. The optical image of the Crab in Eigure[2 
demonstrates a similar discrepancy assuming that the above mentioned optical 
”counter-jet” is the actual NW jet. 

Unfortunately, the compact Vela PWN has not been detected in the optical 
despite considerable efforts (Mignani et al. 2003 Moran et al. [2014) hence 


direct comparison with the Crab PWN is not possible in this band. 


Zynzin[ 


et al. (2013) recently reported the extended feature seen in the NIR (K^ band) 
which could be associated with bar at the base of the SE jet or could be analog 
of the Crab’s inner knot. However, this result still requires confirmation. 

It is also interesting to contrast the radio morphologies of the Crab and 
Vela PWNe. The filamentary structure of the Crab resembles that of Vela X, 
however, the latter has a much larger angular extent (r ~ 60' for Vela vs. 
^ 2.2') and it is much more asymmetric. On the other hand, the Vela PWN is 
a factor of 7 closer (hence, it should appear larger) and a factor of 10 — 20 older 
(hence, it had more time to expand). Therefore, the different angular sizes are 
not surprising. We also note that in the Vela PWN the TeV emission comes 
from the region of brightest radio filament ( [Trail et al.[[1997 Aharonian et al. 


2006b) which is filled with the ejecta (based on the X-ray spectra; LaMassa 
et al.[[20Q8 ) thus providing denser target for the putative relativistic hadrons 


that might be present in the pulsar wind. Therefore, it may turn out that 
some of the prominent thermal Crab filaments are TeV bright. If confirmed, 
it could provide evidence for the elusive hadronic component in pulsar winds 


Arons (1998). Unfortunately, current resolution of the HESS and VERITAS 


telescopes does not allow to test this hypothesis (the existing data only suggest 


that the TeV emission is confined to within <1.7' from the pulsar; Aharonian 


et al.[[2QQM ). Although it is plausible that the filamentary radio morphology 
from the Crab PWN is analogous to that of Vela X, there is no analogy in 


the Crab for the compact radio PWN found in Vela by Dodson et al. (2003). 
The ATCA images reveal a double lobe structure, with the lobes being on 
each side of the X-ray PWN symmetry axis. The radio lobes, extending out 
to 3' — 4', exceed the size of the X-ray arcs by a factor of 5 but nonetheless 


In the Radhakrishnan and Deshpande ( 2Q0l[ > model the arcs are the traces of the particle 
beams from the two magnetic poles and the inner jets are the Doppler-boosted projections 
of the beams. 
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Fig. 6 Photon index maps for Crab ( [Mori et al.|2Q04| > and Vela PWNe ( [Kargaltsev et al.| 
|2013a| obtained from Chandra ACTS data. The color bar at the bottom shows photon index 
values. 


they appear to be filled with the faint X-ray emission which is well seen in 
the harder (1—8 keV) band (see Figure [^. Even if a similar structure in the 
Crab would be smaller by a factor of 10 — 100, it should have been resolved 
in the JVLA images (unless it is buried under the much brighter filamentary 
structure). We also note that no wisp-like structures are seen in the X-ray, 
optical, or radio images of the Vela PWN. 

The X-ray spectral indices of the Crab and Vela PWNe are very different. 
Figure shows spectral maps for the Crab and Vela PWNe. One can see that 
the spectra are the hardest (photon indices are the smallest) for the inner ring 
(in Crab PWN) or for the arcs (in Vela PWN) suggesting that these struc¬ 
tures are associated with freshly injected accelerated particles. Interestingly, 
for the outer arc in Vela PWN the spectrum softens noticeably away from the 
symmetry axis while this is not the case for the inner arc in the Vela PWN or 
the inner ring in the Crab PWN. 
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1.3 Bow shocks and tails: PWNe around supersonically moving pulsars. 


Pulsar average 3D velocities have been found to be ^ 400 km s“^ for isotropic 
velocity distribution (see Hobbs et al.|2005 ). This implies that the majority of 
pulsars stay within their host SNR environment for a few tens of thousands 
years although some particularly fast-moving pulsars can leave it earlier. Once 
the pulsar leaves the SNRp^ it moves in a very different environment which 
has a much lower sound speed. For comparison, the sound speed in the middle- 
aged SNR(such as Vela SNR) can be on the order of a few hundred km s“^. 
The transition between the two very different environments should have a dra¬ 
matic effect on the PWN of the high-speed pulsar. Once the pulsar is moving 
in the medium where its speed substantially exceeds the ambient sound speed 
(i.e. Mach number J\A = v^jcg ^ 1, where and Cq are the pulsar and sound 
speeds, respectively), the PWN shape should be strongly distorted by the ram 
pressure of the medium. If initially the wind was isotropic (this is obviously 
a great oversimplification, see above) the PWN would acquire a cometary 
shape with the pulsar wind being confined to within the surface formed by 
the contact discontinuity (CD) separating the shocked ambient medium and 
the shocked pulsar wind (see Fig. [^. Typically, it is assumed that the pulsar 
wind is shocked in the termination shock (TS) at distances, r^, substantially 
smaller than the distance to the contact discontinuity, Ted, even at the apex 
of the bowshock (see Figure [^. For very fast moving pulsars Ved may become 
so small that for some of the X-ray emitting electrons the gyration radius, r^, 
would become ^ ^ Ved which may lead to leakage of the electrons from the 

apex of the bowshock ( Bandiera||2QQ8 see also below). Numerical simulations 
by Bucciantini et al. ( 2QQ5[ ) indicate high flow speeds Vp) in the shocked 
pulsar wind outflow behind the moving pulsar, suggesting that an extended 
pulsar wind tail should fornj^ For realistically anisotropic pulsar winds (with 
equatorial and polar outflows), in addition to the Mach number, the appear¬ 
ance of the head of the bow shock PWNe and properties of the pulsar tails 
(to within a few from the pulsar) should also depend on the angle between 
the velocity vector and the spin axis of the pulsar. These effects have been in¬ 
vestigated numerically by Vigelius et al. (2007) in the limit of non-relativistic 
3D hydrodynamics who found that the bow shock morphology is only weakly 
affected by the pulsar wind momentum flux anisotropy but the morphology 
of the pulsar wind flow in the tail is strongly affecteef^ On the other hand, 
the ambient medium non-uniformity was found to be greatly affecting the bow 


Alternatively, an o ld SNR can break-up and dissolve. 


Bucciantini et al. 


(20051 do not extend further than 25rcd,o (where 
of the bowshock) from the pulsar due to numerical 


The simulations of 

^cd,o is the scant-off distance at the apex 
challenges. Also, the model neglects the impact of the magnetic field on the flow dynamics. 
It is reasonable to expect that the pulsar tail physics may have some similarities with that 
of leptonic AGN jets for the case when the pulsar spin axis are parallel to its velocity vector 
(except that magnetic hoop stress may turn out to be larger in the case of pulsar tails). 
Therefore, some of the AGN jet simulations may be relevant for the pulsar tails. 

The simulations by Vigelius et al. (20071 only extended for a few rcd,o from the pulsar, 
much smaller than the scales of extended tails seen in X-rays). 
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Fig. 7 A schematic representation of PWNe around the pulsar at rest (left) and superson¬ 
ically moving pulsar (right) for an idealized case of isotropic pulsar wind. 


shock symmetry and shape. Overall, Vigelius et al. (2007) concluded that “the 
anisotropy of the wind momentum flux alone cannot explain the observed bow 
shock morphologies”. The simulations also show that Kelvin-Helmholtz (KH) 
instabilities can develop if the ambient medium exhibits a large pressure gra¬ 
dient. These could be further amplified if the relativistic nature of the pulsar 


wind flow is taken into account due to the increased velocity shear (Bucciantini 
et al.||2005 ). In such situation it is possible that the shocked ambient material 


can be entrained in the pulsar wind flow altering its structure, dynamics, and 
emission properties. The entrainment of ambient matter in the pulsar wind is 
largely an unexplored area (see, however, Lyutikov|2QQ3 Morlino et al.|2Ql^ . 
Further simulations of bow shock PWNe combining these effects (relativis¬ 
tic velocities, ambient medium non-uniformity and entrainment, pulsar wind 
anisotropy and dynamical role of the magnetic fields, 3D geometry and insta¬ 
bilities) can provide a realistic picture for comparison with the observations. 
It may also be possible to make progress by advancing the analytical models 
of these outflows. Romanova et al. (2005) has constructed a model of a pulsar 
magnetotail for the axisymmetric case (the pulsar velocity is co-aligned with 
its spin axis). In this model the pulsar wind remains collimated at large dis¬ 
tances from the pulsar forming a magnetotail where an equipartition is reached 
between the magnetic energy and the relativistic particle energy. The model 
predictions for the shape of the magnetotail appear to agree with the data 
in some cases (e.g., PSR JllOl-6101; Halpern et al. 2014 and PSR J1747- 
2958; [Gaensler et al. 2004 Yusef-Zadeh and Gaensler 2005) and disagree in 
others (e.g., PSR J1509-5850; Kargaltsev et al.|2008 ), possibly, discriminating 
between the axisymmetric and non-axisymmetric cases. 

Chandra and XMM-Newton observations provided for the first time X-ray 
images of PWNe around supersonically moving pulsars (see examples in Fig- 
nres [8| and 11). Some of these images display structures (see e.g., Figs.|^ 
andl^ that are broadly consistent with the theoretical expectations. Indeed, 
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Fig. 8 X-ray (red) and radio (blue) images of pulsar tails produ ced from archiv al Chan¬ 
dra and JVLA data. The J1509-5850 radio image is based on the ( |Ng et al.|2Q10| > analysis 
of Australia Telescope Compact Array observation. 



Fig. 9 Chandra ACIS images of the head regions of B0355+54, J1509-5058, and Mouse 
PWNe (left to right). Notice very different morphologies of B0355+54 and J1509-5058 
PWNe. The images are produced from archival Chandra data. 


in the X-ray images shown in Fig. [^one can identify bright PWN heads ac¬ 
companied by much fainter extended tails (see Fig.|^. Interestingly, the X-ray 
bright PWN head may or may not be bright in the radio [c.f. PWN of PSR 
J1747-2958 (“Mouse”) and PWN of PSR J1509-5058]. Radio polarimetry of 


two extended tails (Mouse and J1509-5058; Yusef-Zadeh and Gaensler 2005 
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and Ng et al.]|2QlQ ) shows that the magnetic field direction is predominantly 
transverse in the case of the J1509-5058 tail and aligned with the tail in the 
case of the Mouse tail. This may suggest that the spin axis is more aligned 
with the velocity vector for J1509-5058 than for the Mouse pulsar (see Fig. 3 
from Romanova et al.|2005 ). Recently, deep, high-resolution observations with 
Chandra revealed interesting structures of the bright bow shock PWN heads. 
Images of J1509-5058 and B0355+54 PWNe, shown in Figure exhibit con¬ 
trasting morphologies. The head of the B0355+54 PWN shows symmetric, 
filled “mushroom cap” morphology with emission being somewhat brighter 
near the center than on the sides. On the other hand, the head of the J1509- 
5058 PWN, comprised of two bent arcs resembling a bow structure (the arcs, 
however, do not quite connect to the pulsar), is mostly “empty” except for the 
slight extension just behind the pulsar. This structure remarkably resembles 
that of the nearby {d ~ 250 pc) Geminga PWN (see Figure [Tq|) as it would be 
seen at a much larger distance of J1509-5058 (d ~ 4 kpc). The bow-shaped 
X-ray emission can either be associated with the forward shock in the ambient 
medium (unlikely, because the pulsar velocity must be very high to produce 
X-rays by heating ISM) or pulsar jets. In the latter scenario the outflows from 
JI509-5058 and Geminga must be dominated by the luminous jets rather than 
the equatorial component (cf. Grab and Vela PWNe). This may be difficult to 
reconcile in the Komissarov and Lyubarsky ( 2004| ) model where the jet forma¬ 
tion is intimately connected to the diverted by the magnetic field hoops stress 
equatorial outflow (backflow) which helps to collimate the polar outflow. Fur¬ 
thermore, the recent 3D simulations (see §2.2 and 2.3) suggest reduced axial 
compression and weaker jets compared to the 2D simulations. On the other 
hand, most numerical simulations are designed to reproduce the Grab and Vela 
structures with a large angle between the NS magnetic dipole and spin axis. 
If this angle is small, the outflow dynamics could be substantially different. 
If the side arcs of JI509 are indeed jets, it would also be difficult to explain 
the ordered helical magnetic field morphology in the extended tail (revealed 
by the radio polarimetry; Ng et al^ 2010) because such a structure would be 
more natural for the axially symmetric case (Romanova et al. 2005). Thus, 
although it is plausible that qualitative morphological differences in the ap¬ 
pearances of compact PWNe can be attributed to the geometrical factors (i.e. 
angles between the line-of-sight, velocity vector, spin axis, and dipole axis), 
these dependencies are yet to be understood. 

Particularly interesting and puzzling is the transition region between the 
bright PWN head and the faint extended tail. For instance, in the B0355+54 
PWN (“Mushroom”) the drop in the surface brightness at the trailing edge of 
the Mushroom “cap” is nearly as sharp as at the leading edge (this makes it 
unlikely to be due to the synchrotron burn-off). In the conventional (isotropic 
pulsar wind) bow-shock tail models this could be associated with the back sur¬ 
face of the bullet-shaped termination shock ( [Bucciantini et ar]|2QQ5[ |Gaensler| 
et al. 2004). However, this interpretation does not appear to work for the 


J1509-5058 PWN which lacks the emission from putative back surface of the 
termination shock. The transition is also much smoother in the Mouse PWN 
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Fig. 10 570-ks Chandra ACIS image of Geminga PWN (Posselt et al. in prep.). 


(see Fig.right panel). The B0355+54 PWN image also shows much fainter, 
narrow “stem” attached to the Mushroom cap which makes it tempting to 
associate the bright trailing edge of the Mushroom cap with the equatorial 
termination shock that has been pushed back by the ram pressure. In this 
scenario the stem and the brighter middle part of the Mushroom cap would 
be associated with polar outflow (a jet). However, even in this case the drop 
in the brightness at the trailing edge of the Mushroom cap may be too abrupt. 
For instance, in the deep Chandra images of the compact Vela PWN one can 
see the effect of the motion onto the inner ring (commonly associated with the 
termination shock) with the particles being blown back off the inner ring (see 
Fig. right panel) . We do not observe such a smooth transition behind the 
Mushroom “cap” (Klingler in prep.). 

The PWNe behind several very fast-moving pulsars display puzzling mor¬ 
phologies (see Figure 11). These are the “Lighthouse nebula” with PSR JllOl- 


6101 ( Halpern et al.pQM Pavan et al.|2Qll Tomsick et al.|2Q12 ), the “Guitar 


nebula” with PSR B2224+65 ( Johnson and Wang|2QlQ Hui and Becker 1 2007 ) 

and the “Turtle nebula” with PSR J0357+3205 ( Marelli et al.|2013 ). The first 
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Fig. 11 Puzzling PWN morphologies are seen in the Chandra images of three high-speed 
pulsars. The arrows show pulsar proper motion directions. In two cases (left panels) the 
extended structures are orthogonal to the pulsar proper motion. For PSR J0357+3205 (right 
panel) the extended X-ray emission brightens further away from the p ulsar while it is very 
dim in the immediate vicinity of the pulsa r (sho wn in the inset). See |Pavan et ah ]( [20l4t ; 
|De Luca et al.| ( |2013| >; pohnson and Wan^ ( |2010] ) for detailed analysis. 


two display bizarre extended features orthogonal to the pulsar’s proper motion 
directions. The third PWN represents a long and luminous tail in the direc¬ 
tion opposite to that of the pulsar’s motion; however, close to the pulsar the 
tail is very faint (undetectable) with no sign of a bright “head” (or compact 
nebula) near the pulsar (cf. Mouse, B0355+54, or J1509-5058 PWNe). Sim¬ 
ilar puzzling behavior is seen for PSR J1101-6101 in the Lighthouse PWN. 
To explain these structures, several hypothesizes have been suggested. Side¬ 
ways structures in PSRs J1101-6103 and B2224+65 could be pulsar jets (e.g., 
Johnson and Wan^|2010[ Pavan et al.||2014 ), although the one-sidedness and 


high X-ray efficiencies {Lx / E) of these structures remain puzzling. The leak¬ 
age of the wind particles from the apex of the bow shock pushed too close to 
the termination shock may be an alternative possibility (Bandiera 2008). In 
the latter case, the morphologies of the sideways features are expected to fol¬ 
low the morphology of the magnetic field in the surrounding ambient medium 
(which, in these cases, is the ISM well outside the pulsar’s host SNRs). 


The bowshock PWNe with tails are often found around older pulsars mov¬ 
ing fast through the rarefied ISM (perhaps due to the observational bias; see 
KP08). Unlike torus-jet PWNe, the bowshocks are often seen in However, 
Ha bowshocks and bright X-ray tails are rarely seen together (exceptions are 
PSR J1741-2054, PSRs J2124-3358, and binary B1957-1-20). The scarcity of 
such cases suggests that the pre-ionization of the oncoming ISM by the high- 
energy radiation from the compact PWN and/or pulsar can play an important 
role. The cases with both He, and radio emission are also very rare but this 
may be due to the limited number of objects observed sufficiently deep in the 
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Fig. 12 ifo; (left) and far-UV (right) images of the bow shock around PSR J0437-4715. 
Faint, amorphous diffuse emission seen around the pulsar in the right panel is the instrumen¬ 
tal artifact (thermal glow of HST/SBC detector) while the structure seen at the bottom is a 
background galaxy. Both right and left panel images show the same area of sky. (Rangelov 
et al. in prep.) 


radio. The sizes of extended pulsar tails can reach 7-8 pc in X-rays (e.g., PWN 
of PSR J1509-5850) and up to ^ 17 pc in radio (e.g., Mouse PWN). These two 
PWNe, well studied both in X-rays and radio, show remarkably contrasting 
behavior. The Mouse PWN appears to be the brightest closer to the pulsar 
both in radio and X-rays while the radio emission from the J1509-5850 PWN 
is lacking near the pulsar and becomes brightest only a few arc minutes away 
from the pulsar (see Fig. [^. This could be attributed to the fact that the 
Mouse PWN is seen more of a head-on but the very long radio tail argues 
against the very small angle between the pulsar velocity and the line-of-sight 
(see Figure 1 in Yusef-Zadeh and Gaensler|2QQ5 ). The morphology of the radio 
PWN of PSR J1101-6101 is similar to that of J1509-5850 radio PWN. Fur¬ 
thermore, the radio polarization measurements indicate that in the Mouse tail 


the magnetic field direction is predominantly parallel to the tail (Yusef-Zadeh 


and Gaensler 

2005 

), in i 

the tail (Ng et al. 

2010) 


2010). The differences may be related to different angles 


between the pulsar velocity and spin axes in these two PWNe. Alternatively 
a different ambient density and entrainment rate could play a role. 

PWNe with bow shocks are particularly interesting because the 
emission allows one to map the structure of the forward shock not only in 
coordinate space but also in velocity space through the measurements of the 
Doppler shifts in hydrogen lines across the forward shock ( Romani et al.|2QlQ 


Brownsberger and Romani 2014). For instance, Romani et al. (2010) performed 


spectroscopic observations for J1741-2054 and measured the radial velocities 
up to ^ 50 km s“^ consistent with the bowshock model implying pul sar speed 
of ^ 150-200 km s“^ and inclination angl^^of about 75° (see also Auchettl| 


17 


The angle between the line of sight and the pulsar velocity vector. 
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Measurements of spectral line {Ha or/and Hi^) fluxes can provide accu¬ 
rate diagnostics of the ambient medium density if other parameters (pulsar 
velocity h^psr, distance to the pulsar d, and the stand-off distance Vcd) are 
well constrained. A nearby millisecond PSR J0437-4715 with E = 5.5 x 
1033 g-i represents an example of such kind. The pulsar, located at 
the (parallax) distance of d = 156 pc and moving with Vpsr = 134 km s“^ 
(for the inclination i = 53° inferred from the bowshock shape), shows a 
prominent Hq, bow shock (Fig. 12) with a stand-off distance = 2.2 x 10^^ 
cm. This implies ISM density uh = E/ATrcrgiTiHV^^^ = 0.1 cm“^. For slow 
(l^sr < 10^ km s“^) bowshocks the yield (i.e., number of 11^^ p hotons 
per incoming 


neut ral) is = 0.6(u/100 kms ^) according to |Heng 


and McCray (2007) and the Ha flux ~ 0.0074(u/134 kms ^)^/^(rg /2.2 x 


10^^ cm)^(d/156 pc) ‘^{uh/^H cm ^)^hi photons cm ^ s ^ (see e.g., Chat- 
terjee and Cordes]|2002| Brownsberger and Romani]|2014 ). The measured H, 


-2 


flux from the J0437-4715’s bow shock apex, = 6.7 x lO”"^ photons cm 
s“^ (Brownsberger and Romani 2014), implies that the neutral H fraction 
~ 0.9. This is somewhat surprising because the classical Stromgren ra¬ 
dius for the ionizing radiation produced by the pulsar appears to be much 
larger than the distance to the Hq, bowshock apex indicating that the clas¬ 
sical formula is probably inapplicable when ionizing radiation is X-rays and 


the ionizing source is moving fast van Kerkwijk and Kulkarni (2001); Lyu- 


tikov et al. (2015). PSR J0437-4715 bow shock is the only one from which the 


far-UV (FUV; see Fig. 12) emission has been detected (other pulsars with Hq, 
bowshocks may simply be too far, so FUV photons are easily absorbed). The 
measurement of the FUV spectral slope, a ^ 1.5 (for E^, oc z^“^), suggests ei¬ 
ther a non-thermal emission mechanism (synchrotron?) with radiation coming 
from the vicinity of the forward shock (FUV and H-alpha bow-shocks coincide 
within the measurements errors, 0.5") or line emission (with multiple un¬ 
resolved lines resulting in an effectively flat PL spectrum). No firm detection 
of X-ray emission has been reported for the pulsar wind of J0437-4715 yet, 
although an analysis of archival Chandra ACIS data indicates that a compact, 
^ ^ 2" — 3", asymmetric (most of the emission is ahead of the pulsar) PWN 
with the luminosity of 3.8 x 10^^ erg s“^ may be present (Rangelov et al. in 
prep.). 


So far B1951-1-32, J1509-5850, and J2124-3358 are the only three solitary 
pulsars where both X-rays and Hq, have been detected. However, in J1509- 
5850, Ha emission is believed to come from the photo-ionized medium rather 
than from the bow shock (see the discussion in Brownsberger and Romani] 


2014) and for B1951-1-32 the surrounding environment is complex due to pos¬ 
sible contribution to Ha from the filaments of the host SNR CTB 80. Therefore, 
it is important to perform deep H-alpha imaging and spectroscopy of nearby 
pulsars with X-ray bow shocks to study the connection between the pulsar 
wind and shocked ambient medium regions. 
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1.4 PWNe in binaries. 


Sufficiently energetic non-accreting binary pulsars may be able to power de¬ 
tectable PWNe. In binary systems, in addition to the pulsar parameters, the 
PWN morphology and appearance would depend on the properties of the bi¬ 
nary companion and parameters of the binary orbit. Even if the companion 
star is lacking any wind itself (e.g., it is a cold white dwarf or old low-B neu¬ 
tron star) the pulsar wind will be affected by the ram pressure due to its own 
orbital motion if the ram pressure caused by the orbital motion is comparable 
to or larger than the ambient pressure (or ram pressure due to the motion 
of the binary as a whole). Thus, binary motion can strongly affect the pulsar 
wind flow (and morphology of the magnetic field) downstream of the termina¬ 


tion shock or even the termination shock itself (see e.g., Bosch-Ramon et al. 
2012 Zabalza et al.||2Ql'3 ). If the companion star has a powerful wind the in¬ 


teraction becomes even more complex with the outcome critically dependent 
on the ratio of momentum fluxes of the two winds rj = Ec ~^/(where 
Mw and are the mass loss rate and the massive star wind velocity), a key 
parameter in colliding wind binaries (see Dubus (2013)). If 77 ^ 1, the pul¬ 
sar wind dominates (e.g., for B1957+20 - the original “black widow” system) 
while in the opposite case the companion star wind will be dynamically dom¬ 
inant. In general, 77 can vary with the orbital position (1) if the pulsar wind 
(or the massive star wind) is anisotropic, and the spin axis of the pulsar (or 
the massive star) do not coincide with the orbital angular momentum vector 
or ( 2 ) if the orbit is highly eccentric. 

The most famous example of a pulsar binary system where all these effects 
play a significant role is LS 2883 with the young energetic pulsar B1259- 
63 (B1259 hereafter) in an eccentric 3.4-year orbit around a massive 0-star. 
Although direct observations of pulsar wind in B1259 may not be feasible 
(except perhaps for the VLBI imaging observations; see Moldon et al.||2Q12 ), 
there are indirect ways to learn about the pulsar and stellar wind proper¬ 
ties and their interaction. These include multi-wavelength spectrum and flux 
measurements as a function of orbital phase (Chernyakova et al. 2014), pul¬ 
sar radio signal variability measurements ( Shannon et al.||2QT4 and references 
therein), and, a high-resolution X-ray imaging, which recently revealed a dy- 

Pavlo^ 


2014b 


namic structure associated with the binary (Kargaltsev et al 
eFahjlM^ . An X-ray emitting cloud was found to be moving away from the 
binary with the velocity of ~ 0.07c ( Pavlov et al.||2015 ), which, together with 


the lack of deceleration, implies either a hadronic cloud with very large mass 
^ 10^^n(d/2.3 kpc)^ g moving in the 0 -star wind with density n or a lighter 
cloud moving in the rapidly expanding, unshocked relativistic pulsar wind (im¬ 
plies 77 ^ 1 in the polar 0-star wind). The former scenario implies that the 
cloud was ejected from the binary during the 2011 periastron passage, when 
the pulsar interacted with the excretion disk of the massive 0-star. However, 
the corresponding kinetic energy of the cloud must be very large, ^ 2 x 1045 
erg, and it must have been launched via a complex interaction between the 
pulsar wind and excretion disk of the 0 -star with the energy source being 
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Fig. 13 Chandra observations (on MJD 55912, MJD 56431, and MJD 56696 left to right) of 
LS 2883/B1259-63 reviling the relativistic, v = (0.074zb0.006)c, motion of cloud ejected from 
the binary. See the corresponding movie at http://home.gwu.edu/~kargaltsev/B1259.html 


problematic (pulsar’s E can only provide ^ 10^^ erg during the disk passage). 
The latter scenario is at odds with the common assumption of ^ 1 for such 
kind of binaries but it does not require extreme values of mass and energy for 
the cloud ( Pavlov et al.||2Q15 ). 

We note that there are other systems where pulsars (albeit these may not 
be as young and energetic) might orbit massive (often Be or O type) stars, 
and it is plausible that other TeV gamma-ray binaries (e.g., LS 5039, LS I -h61 
303, and HESS J0632+057) also host pulsars (|Dubus||2013 Kargaltsev et al 


2014a and references therein). Indeed, Durant et al.| ( |20TT ) reported evidence 


for amorphous arcminute-scale X-ray emission with a hard spectrum around 
LS 5039 which was interpreted as synchrotron radiation from ultrarelativistic 
(pulsar wind?) particles escaping from the system. 

Among other types of binaries which can shed light on the properties of 
pulsar winds at smaller distances from the pulsar through the interaction with 
the companion are the famous double pulsar (see e.g., Arons et al. 2005), 


recently reported very eccentric binary with PSR J2032+4127 (Lyne et al. 
2015), and black widow pulsar B1957+20 ( Huang et al.||2012 ). 


1.5 PWN spectra 


It is commonly accepted that radio-to-MeV emission in PWNe is the syn¬ 
chrotron emission. This naturally explains the observed high degree of po¬ 
larization in X-rays, optical, and radio (Moran et al. 2013| Bietenholz et al.| 
2004 Weisskopf et al.||1978 ). Multiwavelength observations and modeling show 


that IC scattering on Cosmic Microwave Background and NIR/IR background 
starts to dominate PWN spectra in GeV (Acero et al. 2013[) and completely 
dominates synchrotron at TeV energies ( Atoyan and Aharonian|p~996[ |Volpi| 
et al.|[2QQ8l [Gelfand et al. 2009[ [Torres et al.||2013[ and |01mi et al.||2014[ s^ 
the Grab and Vela PWN spectra in Fig. 14 for example). The theoretical mod¬ 
els of particle acceleration and non-thermal emission in PWNe predict, in the 
ideal MHD framework, a cutoff around 150 MeV (see §3.1). These models well 
reproduce the observed cutoff in the Grab Nebula spectrum (see left side of 
Fig. 14). However, the slow acceleration rates of ideal MHD models (as, e.g.. 
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20081 and Vela PWN 


Fig. 14 Multiwavelength spectra of the Crab (left; from 
(right; spectrum from r = 6' from pulsar shown in color while the relic Vela X plerion 
spectrum is in grey; from |Mattana et al.|20lH . 


Table 1 Photon indices (P) and electron SED slopes (p = 2P — 1) measured from Chandra 
ACIS data for the innermost regions in 9 bright PWNe ( [Kargaltsev et al.|2014c| >. 


PWN 

r 

P 

Crab 

1.80 ±0.05 

2.6 

Vela 

1.30 ±0.05 

1.6 

3C58 

1.9 ±0.07 

2.8 

G320.4-1.2 

1.4 ±0.1 

1.8 

Kes 75 

1.9 ±0.1 

2.8 

G21.5-0.9 

1.40±0.06 

1.8 

Gll.2-0.3 

1.5 ±0.1 

2.0 

CTB 80 

1.7 ±0.1 

2.4 

G54.1+0.3 

1.50 ±0.05 

2.0 


diffusive shock acceleration) do not allow fast variability in the nebular emis¬ 
sion (see §3 for Crab flares). The questions about the location and distribution 
of acceleration sites as well as acceleration mechanism in the pulsar winds and 
PWNe also remain open (see §2.4). 


If the pulsar wind contains relativistic protons, it is possible that hadronic 
emission (due to neutral pion decay) contribution can become appreciable in 
dense environments ( Bartko and Bednarek|[^Q8| ). Multiwavelength emission 
from bow shock PWNe (including pulsar tails) produced by pulsars moving 
in a low-density ISM (outside their host SNRs) should be purely leptonic. 
However, surprisingly few of these objects have been detected in TeV (one of 
the deepest limits, 0.1% of ^ in 1-10 TeV, is obtained with VERITAS for the 
tail of B0355+54; Brett McArthur, private communication). Recent review of 
the observational X-ray and TeV properties for the population of 91 PWNe 
can be found in Kargaltsev et al. (2013b). 
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1.5.1 Spatially resolved X-ray speetra 

Deep Chandra ACIS observations of a few bright PWNe allow us to create 
spectral maps with high-spatial resolution ( Kargaltsev et aL]|2Q13a[ ). These 
maps are expected to provide manifestation of synchrotron burn-off (X-ray 
spectra should become softer farther away from the pulsar) which depends 
on the strength of magnetic field and bulk flow speed (Kennel and Coroniti 
1984b Wang et al. 2001 Reynolds [2009 ) but they may also contain signa¬ 


tures of spatially distributed (in-situ) particle acceleration or rapid particle 
diffusion. These spectral maps demonstrate that the pulsar spectra measured 
just downstream of the termination shock can differ substantially. From Table 
[^one can see that although for the Crab PWN the inferred (assuming syn¬ 
chrotron emission model) slope of the electron SED, p = 2T — 1, is consistent 
with the p = 2.1 — 2.8 expected from the commonly invoked Fermi accelera¬ 
tion mechanism (Achterberg et al. |2QQ1[ |Sironi and Spitkovsky|2QQ9a [|2Qllb[), 


for the Vela PWN the spectrum is much harder {p ^ 1.6) suggesting that 


a different mechanism might be at work (e.g., magnetic reconnection; Sironi 
and Spitkovsky||2Q14[ ). At least some of the long pulsar tails (e.g., tail of PSR 


J1509-5850) tend to show very little evidence of cooling (in terms of spectral 
softening in X-rays) which either suggests ongoing in-situ acceleration along 
the tail or extremely fast bulk flow (Klingler et al. in prep.). 


1.5.2 X-ray and TeV effieieneies of PWNe 


The substantial number of PWNe detected in X-rays and TeV allows one to 
investigate the population properties. Here we will only consider X-ra y and 
TeV radiative PWN efficiencies and refer the readers to 


Kar¬ 


galtsev et al. (2013b) for the analysis of other TeV and X-ray properties of 
PWNe. Figure |1 5 (based on the information collected in Tables 1-3 in 


Kar- 


galtsev et al.|2Q13b{ with some updates) shows the X-ray and TeV luminosities 


of PWNe (and PWN candidates for TeV). Notice a very large spread of X- 
ray efficiencies and a noticeably smaller spread for the TeV efficiencies. While 
there is a clear correlation between E and px there is no noticeable correla¬ 
tion in TeV. Finally, the majority of PWNe that are underluminous in X-rays 
appears to be around y-ray loud, radio-emitting pulsars with small magnetic 
inclination angles (based on Rookyard et al.|[2Q14 ). 


1.5.3 Seareh for new PWNe in GeV j-rays. 

PWNe are the most numerous source class that emerged from the H.E.S.S. 
Galactic Plane Survey ( Carrigan et al.|2Q13 ), and about 100 PWNe or PWNe 
candidates are known in X-ray. Nevertheless only very few sources are detected 
in the MeV-GeV range, because this energy range falls between the tail of 
the synchrotron emission and the rising part of the IG emission and because 
of the lack of sensitive MeV instrumentation with good angular resolution. 
Often, search for PWN at GeV energies has to be carried out in the presence 


































Pulsar-Wind Nebulae 


25 



33 34 35 36 37 38 39 

Log Edot 



Fig. 15 Top: X-ray luminosities of PWNe and PWN candidates vs. pulsar’s E. TeV PWNe 
and TeV PWN candidates are shown in red. The dotted straight lines correspond to constant 
X-ray efficiencies; the u pper bound, logL^^ = 1.51 log ^ — 21.4, is shown by a dashed line 
( [Kargaltsev et alT| 2012 | ). The PWNe detected in GeV by Fermi are marked by stars. Blue 
filled circles are the pulsars with confidently measured small (< 10 °) magnetic inclination 
angles from |Rookyard et al.| ( [20T4| . Bottom: TeV luminosities of PWNe and PWN candidates 
vs. pulsar’s E. Thin error bars mark questionable associations. The PWNe undetected in 
X-rays are marked by circles. PWNe detected by Fermi are marked by stars. The dotted 
lines correspond to constant values of the TeV 7 -ray efficiency rjj = L^jE. The detection 
of TeV emission from Geminga region is, so far, based solely on Milagro result [Abdo et aT] 
( |20Q9| > which has not been confirmed by any other observatories (e.g., VERITAS). 
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of 7 -ray-loud pulsar by looking for the off-pulse emission which could come 
from the PWN (see jAbdo et al.|2QT3 ). In Acero et al. (2013), a search for GeV 
emission from 58 TeV PWNe and unidentified sources was performed, with the 
requirements of (1) good GeV and TeV spectrum connection and (2) extended 
emission. A total of thirty sources were detected in the GeV range, for energies 
above 10 GeV; among them, 11 sources are PWN candidates, and three are 
reliably identified as new PWNe. These new sources are associated with young 
(age between 1 and 30 kyr) and powerful pulsars with E between 10^^ and 
1039 g-i interesting to study how the multi wavelength properties 

of the PWNe evolve with the properties of the host pulsar. It was found in 


Acero et al. (2013) that there is no correlation between the GeV luminosity 


and the age and the spin-down luminosity of the pulsars, and the same for the 
GeV-to-TeV luminosity ratio. On the other hand, the GeV-to-X-ray luminosity 
ratio appears to increase with age in agreement with theoretical models (see 
e.g., Torres et al.||2013 ). Even less is known about PWN emission at the MeV 
energies due to the poor imaging capabilities of the MeV telescopes. The dip 
is expected to occur at these energies in the spectra of ^10 to ^100 kyr-old 
PWNe which are most frequently found in X-rays and TeV. Measuring the cut¬ 
off energy of the synchrotron spectrum in hard X-rays constrains maximum 
energies of the electrons in PWNe. 


1 . 5.4 Multiwavelength speetra 

Spatially-resolved spectroscopy is primarily done in X-rays, thanks to the su¬ 
perb resolution of Chandra. It remains rather challenging to obtain spatially- 
resolved spectral measurements at other frequencies. In the optical-NIR this 
is primarily due to the faintness of the PWN emission (except for the Grab 
PWN) and contamination by various background sources. In the radio, this is 
challenging because of the interferometric nature of high-resolution observa¬ 
tions and difficulties in combing the requirements of imaging both the large- 
and small-scale structures. At higher energies (MeV, GeV and TeV), the res¬ 
olution of the existing telescopes is often insufficient to resolve PWNe, and in 
those cases when large relic PWNe can be resolved ( [Aharonian et al.| 2006e 


Abramowski et al. 2012) the limited signal-to-noise ratio typically precludes 


spatially resolved studies (see, however, Aharonian et al. 2006a Van Etten 


and Romani||20ir ). Therefore, the MW spectra of most PWNe are, by neces¬ 


sity, spatially-integrated which in many cases introduces bias and systematic 
uncertainties that are difficult to account for ( Van Etten and Romani|[2011 ). 
Indeed, as one can see from Eig. the MW spectrum of the compact Vela 
PWN is very different from the MW spectrum of the relic plerion Vela X 
(which may be a peculiar one in other respects also, see below). The PWN 
luminosities shown in Eig. are calculated for very different regions in most 
PWNe and therefore should be treated as such (i.e. cannot be treated as if 
both X-ray and TeV photons are emitted from the same PWN region). If both 
TeV and X-ray efficiencies (luminosities) of PWNe are compared to those de¬ 
rived from a one-zone PWN model, one should be aware of the limitations 
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of this approach. Multi-zone models, taking into account both advection and 
diffusion, appear to be the next logical step in theoretical development and in 


preparation for Cherenkov Array Telescope (CTA; see e.g., de Oha-Wilhelmi 
et al.|[2M3 Fortson and CTA Consortinm||2Q15 ). 


As an example (albeit perhaps an unusual one) of a relic PWN one can 
consider Vela X. The nearby Vela SNR 8° in diameter) has a large region 


of non-thermal radio emission surrounding the Vela pulsar (see e.g., Bock 
et al.lp^oM ). One of the brightest radio filaments in Vela X, positioned at the 


southwest of the pul sar, was detected in X rays with ROSAT, ASCA, Suzaku. 


and XMM-Newton (Markwardt and Ogelman 


1995 


1997 


Mori et al. 


2008 


LaMassa et al. 2008) and, more recently, in GeV y-rays with AGILE ( [Pellizzoni 


et al.|201Q ) and Fermi LAT ( Abdo et al.|2010b ) and very high energies (0.5-70 


TeV) with HESS ( Aharonian et al.||2006c ) and CANGAROO (Enomoto et al 


2006| ). The bright X-ray and VHE emission regions are positionally coincident 


(they sometimes referred to as a “cocoon”), and have been commonly dubbed 
a relic PWN displaced to the south by the unequal pressure of the reverse 
shock propagating within the SNR. However, subsequent deeper observations 
with H.E.S.S. ( Abramowski et al.|2012 ) and Fermi LAT ( Grondin et al.|20f3 ) 
revealed a fainter extended emission whose morphology appears to correlate 
with the the double-lobe, large (^ 2° in extent) structure found at 61 GHz 
in the WMAP images. This yields strong support to the scenario where two 
different populations of electrons are needed to reproduce the radio/GeV halo 
and the X-ray/TeV cocoon, respectively ( de Jager et al.||2QQ8 ). 


2 Theory of PWNe 


2.1 


General properties of pulsar winds and structure of PWNe 


Theoretical studies of PWNe concentrated on interaction of the relativistic 
pulsar wind with the surrounding plasma. The morphology of nebulae is de¬ 
scribed in the scope of magnetohydro dynamic (MHD) models, which proved 
to be very successful. In order to understand the radiation spectra, one has 
to develop realistic particle acceleration models; this is a much more difficult 
task. Here we shortly review the recent development of the PWNe theory; for 
more comprehensive reviews, see, e.g. Kirk et al. (2009); Arons (2012); Buhlerj 
and Blandford (2014). 


The general features of PWNe are basically dictated by the physics of the 
pulsar wind: 

— The pulsar wind is composed predominantly of electron-positron pairs, may 
be with some admixture of ions. The pair content of PWNe suggests that 
the pair density in the wind is rather high, more than enough to ensure 


that the wind could be described as an MHD flow ( |de Jager et al.||1996 
Bucciantini et al^|2011 ). 


— The wind is strongly magnetized; most of the energy is transferred, at least 
initially, as Poynting flux. 
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— The wind is highly anisotropic: the Poynting flux is maximal in the equa¬ 


torial belt and goes to zero at the axis ( [Michel 
Tchekhovskoy et 'ar]|2Q13 ). 


1973 Bogovalov 1999 


— An obliquely rotating magnetosphere produces variable electromagnetic 
fields that propagate in the wind as MHD waves; specifically in the equa¬ 
torial belt, where the magnetic field changes sign every half of period, the 
so called striped wind is formed. 

The PWN is in fact a bubble filled predominantly by relativistic particles 
and magnetic fields; it is inflated by the pulsar wind that continuously pumps 
into the surrounding medium the energy in the form of relativistic particles 
and magnetic fields. Within the nebula, the fields and the particles are roughly 
in equipartition therefore the main question is how the Poynting flux in the 
pulsar wind is converted into the energy of particles (the so called cr-problem). 

Even though the details of the energy transformation process remain ob¬ 
scure, the general picture is robust: most of the energy is transferred in the 
wind by alternating electro-magnetic field; and the energy is transferred to the 
plasma when the alternating field dissipates. This conclusion follows from the 
strong anisotropy of the MHD wind, in which energy is predominantly trans¬ 
ferred in the equatorial belt where the striped wind is formed. Of course the 
fraction of the energy transferred by alternating fields depends on the angle 
between the magnetic and rotational axes of the pulsar: an aligned rotator does 
not produce alternating fields at all whereas the energy from a perpendicular 
rotator is totally transferred by alternating fields. Due to the strong anisotropy 
of the pulsar wind, most of the energy is transferred by alternating fields even 


in a moderately oblique case. Komissarov (2013) calculated the fraction of the 


energy flux due to alternating fields making use of the split monopole model 
of the pulsar wind ( Michel||1973 BogovalQ^|1999 ), in which the Poynting flux 
is distributed as sin^ where 0 is the polar angle. He found that even if the 
angle between the rotational and magnetic axes is 45^, as much as 72% of 
the total energy is transferred by alternating fields. In real pulsar winds, this 
fraction is even larger because the Poynting flux in the wind from a rotating 
dipole magnetosphere is more concentrated to the equatorial belt than in the 
monop ole wind; according to numerical simulations by [Tchekhovskoy et al 


i, 


(2013k, the angular distribution of the Poynting flux is close to sin 0 in this 


case. 

Even though the dissipation mechanisms for variable fields in pulsar winds 
are still debated, one has to stress that the particle Larmor radii within the 
nebula significantly exceed the wavelength of the waves in the pulsar wind. 
Therefore these waves could not penetrate into the nebula in any case; even 
if they survive within the wind, they dissipate at the termination shock front 


( Lyubarsk^[2QQ3b Petri and Lyubarsky[[2QQ7 Sironi and Spitkovsky[[2QlTa ). 
Therefore the post-shock flow must be the same as it would be if the dissipation 
has already been fully completed in the wind. The structure of the nebula is 
determined by the distribution of the total energy and mean magnetic field in 
the wind. In the equatorial belt, where most of the energy is transferred, the 
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mean field is weak therefore relatively weakly magnetized plasma is injected 
into the nebula in this region. At higher latitudes, the magnetic field does not 
change sign and the variable magnetic fields could propagate in the form of 
fast magnetosonic waves. These waves efficiently decay as a result of non-linear 
steepening ( Lyubarsky|2QQ3a ) however, they could not transfer a large fraction 
of the Poynting fiux so that the flow magnetization in this region remains large 
even after they decay. 

Taking into account the above properties of the pulsar wind, the ob¬ 
served morphology of PWNe, and first of all the characteristic disk-jet struc¬ 
ture ( Weisskopf et al.|[2QQQ Gaensler and Slan^|2QQ6 ), is naturally explained 
(Lyubarsky 2002). Namely, the disk is formed by the relatively weakly mag¬ 


netized equatorial flow, which transfers most of the energy. In the Crab, such 
a disk is clearly seen in the Chandra map because the X-ray emitting elec¬ 
trons rapidly loose their energy thus making X-rays a tracer for the freshly 
injected plasma. In other PWNe, like Vela, the disk may even not be seen at 
all because the high energy electrons All a much larger volume. The high lati¬ 
tude flow remains magnetized therefore it is compressed by the magnetic hoop 
stress to form a jet-like feature at the axis. An important point is that in the 
highly relativistic, super-sonic (more exactly, super-fast-magnetosonic) wind, 
the magnetic collimation is inefficient; the ”jet” is formed beyond the termi¬ 
nation shock where the flow is decelerated. Inasmuch as the pulsar wind is 
anisotropic, the termination shock is highly non-spherical: it lies much closer 
to the pulsar in the polar regions than in the equatorial belt; therefore the 
”jet” appears to originate from the pulsar (see also Bogovalov and Khangou- 


lian (2002)) 


This general conjecture has been confirmed by axisymmetric MHD simu¬ 
lations ( Komissarov and Lyubarsky||2004 Del Zanna et al. 2004 2006 Volpi 
et al.|20Q8 Bucciantini et al. 2011). These simulations were able to explain also 
nontrivial features of the fine structure, such as a mysterious knot, which is 
located within 1” from the Crab pulsar ( Hester et al.|199^ ). Namely, the knot 
is a Doppler-beamed emission from the patch of the highly oblique termina¬ 
tion shock where the post-shock flow is still relativistic and directed towards 
the observer (Komissarov and Lyubarsky 2004). Simulations with a better 


resolution and therefore with a lower numerical viscosity ( Camus et al.|[2QQ9 ) 
revealed bright fine filaments moving away from the termination shock with a 
good fraction of the speed of light. These are highly reminiscent of the so-called 
wisps of Crab nebula ( Hester et al.||1995 ). 


In spite of all these successes, important basic problems have remained un¬ 
resolved. Namely, axisymmetric simulations reproduce the observed structure 
of the nebula only if the wind magnetization at high latitudes was chosen to 
be relatively small, significantly smaller than one could expect from theoreti¬ 
cal considerations. This discrepancy was resolved only recently when fully 3D 


simulations have been performed (Forth et al. 2014). Let us describe these 
recent developments in more details. 
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2.2 The cr-problem and 3D simulations of PWNe 


The problem of the magnetic to the plasma energy transformation in pulsar 
winds is generally referred to as the cr-problem because the flow magnetization, 
defined as the ratio of the Poynting to the plasma energy fluxes, is typically 
denoted by cr. The pulsar wind starts as a highly Poynting dominated (cr ^ 
10^-10^); on the other hand, there is a pervasive belief that one can account 
for the morphology of PWNe, including the remarkable jet-torus structure, 
only if just upstream of the termination shock, cr does not exceed 0.01. Such 
a tremendous drop in the flow magnetization looks so mysterious that the 
problem was sometimes referred to as the cr-paradox. However, one has to 
stress that what we really need to consider is the mean field because alternating 
fields inevitably decay. As it was mentioned above, they transfer most of the 
energy in the pulsar wind therefore the magnetization due to the mean field is 
not large, which makes the cr-problem not so severe. Let us discuss the issue 
in a bit more details. 


First of all one has to stress that the strong constraints on the wind mag¬ 
netization at the termination shock mentioned above were obtained in spheri- 


cally or axially symmetric models of PWNs (Rees and Gunn|1974 

Kennel and 

Coroniti|| 1984a Emmering and Chevalier 

1987 

Begelman and Li 

1992). The 


reason for the required low value of cr is that in these models, the magnetic 
field strength grows with radius within the nebula so that the field would ex¬ 
ceed the equipartition value and pinch the flow too much if the magnetization 
at the termination shock is not extremely small. The behavior of the mag¬ 
netic field could be easily understood if one takes into account that the field 
in the far zone of the wind is practically azimuthal, and in the axisymmetric 
flow, the field lines remain coaxial circular loops. The radius of the field line 
increases when the flow expands, the field strength being determined by the 
conservation of the magnetic flux within the toroidal magnetic tube. At the 
termination shock, the flow compresses so that the magnetic field increases 
three times. The flow within the nebula is subsonic therefore the pressure and 
the density of the plasma do not change significantly. Therefore the volume of 
the toroidal magnetic tube remains roughly constant. In this case, the cross 
section of the tube decreases when the tube radius increases, which implies an 
increase in the magnetic field roughly linearly with the radius. 


Taking into account that the size of the nebula is about an order of mag¬ 
nitude larger than the radius of the termination shock, one finds that in the 
axisymmetric flow, the field strength in the main body of the nebula exceeds 
that in the wind just upstream of the shock ^ 3 x 10 = 30 times whereas a 
grows roughly three orders of magnitude. The problem can be alleviated if the 


kink instability destroys the concentric field structure in the nebula (Begelman 


1998). Then the magnetic loops could come apart and one expects that in 3D, 


the mean field strength is not amplified much by expansion of the flow. In 
this case, a just upstream of the termination shock might not need to be so 
unreasonably small as was found in axisymmetric models. 
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This idea can be checked only by 3D simulations of plasma flow within 
the nebula. As the first step, Mizuno et al. (2011) simulated the 3D evolution 


of a simple cylindrical model of PWNe developed earlier by [Begelman and| 


Li 


(1992). This model describes a static cylindrical configuration with a rela- 


tivistically hot plasma such that the thermal pressure is balanced by the hoop 
stress of a purely toroidal magnetic field. The simulations clearly demonstrated 
that the kink instability does develop in the system and destroys the regular 
concentric structure of the magnetic field thus relaxing the hoop stress and 
triggering magnetic dissipation. This proves that 3D effects play crucial role 
in the evolution of PWNe. However, these simulations do not claim to model 
PWN, simply because the continuous injection of magnetic flux and energy 
into PWN by their pulsar winds is not accounted for. 

The first realistic 3D simulations of PWNe were performed by [Port h et al.| 


(2014). They used qualitatively the same setup as in 2D simulations, namely. 


the nebula is pumped by a strongly anisotropic pulsar wind with the mag¬ 
netization determined only by the mean magnetic field as if the alternating 
component of the field has completely dissipated. The difference was in the 
magnetization at high latitudes, where the magnetic field does not change 
sign. In 2D simulations, the observed morphology was reproduced only if the 
high latitude a was chosen to be as small as 0.1 even though according to the 
pulsar wind theory, it should remain significant, not less than a few. |Por th\ 
et al. (2014) took a = 1 ^ 3 at high latitudes, as it should be. 


According to the results of their 3D simulations, the azimuthal component 
of the magnetic field is still dominant in the inner part of the nebula, which 
is filled mainly with freshly injected plasma. The hoop stress of this field is 
still capable of producing noticeable axial compression close to the termina¬ 
tion shock and driving polar outflows, required to explain the Crab jet, and 
jets of other PWNe (Fig. 16). However, these are much more moderate than 
in 2D models. In the main body of the nebula, the highly organized coax¬ 
ial configuration of magnetic field is largely destroyed by the kink instability 
(Fig.[^ therefore the global evolution of the PWN in 2D and 3D cases differs 
radically (Fig. [^. If the high latitude magnetization is large, the 2D models 
develop extremely strong polar jets, which burst through the supernova shell. 
In contrast, in the 3D models the ^-pinch configuration is destroyed by the 
kink instability so that the polar outflows are less powerful and eventually lose 
collimation, as observed. 

Important observational constraints are imposed by polarization measure¬ 
ments that reveal a high degree of polarization in the central part of the Crab 
nebula with the polarization vector parallel to the pulsar axis, which may 
be considered as direct evidence for the azimuthal field (Hester 2008). The 


simulations do show that in spite of the strong disruption of the azimuthal 
magnetic field, the polarization remains substantial, particularly in the inner 
part of the nebula. The polarization direction on the scale of the torus clearly 
indicates an azimuthal field because the photon magnetic vector appears to 
curve around the torus. 
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Forth et al. 


(20141. 


Fig. 16 Formation of the polar outflow in the 3D MHD simulations by 
The colour images show the distribution of Ig /3 in the yz plane, where /5 is the ratio of the 
thermal to the magnetic pressure. The black lines show the momentary streamlines and the 
red line the termination shock. 


The 3D simulations also show a bright knot, which was discovered in the 
previous 2D simulations ( Komissarov and Lynbarsky|2QQ4 ) and identified with 
the inner knot observed in the Crab nebula (Hester et al. 1995). This emis¬ 
sion comes from the immediate vicinity of the termination shock, it is highly 
Doppler beamed and originates in the high-speed part of the post-shock flow. 
A correlation was found between the knot position and the flux, such that 
brighter states correspond to a smaller offset between the knot and the lo¬ 
cation of the pulsar, which is in excellent agreement with the recent optical 
observations of Crab’s knot (Moran et al. 2013). The simulated polarization 
degree and polarization angle in the knot also agree with observations. 
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Fig. 17 Field lines in the 3D simulations by 


(2014). The lines are coloured 


according to their orientation, sections with dominating azimuthal component being blue 
and those with dominating poloidal component being red. The surface of the termination 
shock is also shown, using the magenta contour. 


The termination shock is found to be unsteady due to an intricate feed¬ 
back mechanism between the shock and the nebula flow. The inhomogeneities, 
formed in the post-shock flow as a result of this variability, appear as wisps 
emitted from the shock location, in a qualitative agreement with the observa¬ 
tions of the Crab nebula. 


In the main body of the nebula, the kink instability not only destroys 
the regular magnetic field structure; the excited turbulence yields efficient 
magnetic dissipation. In simulations, this occurs at the grid scale via numer¬ 
ical resistivity. To become efficient, it requires creation of ever smaller scale 
structures in the magnetic field distribution. However, it is important that 
the processes which drive the development of such structures occurs on scales 
above the grid scale so that the dissipation rate is determined sufficiently ac¬ 
curately. The simulations with doubled resolution show the same dissipation 
rate, which suggests that the high degree of dissipation observed in these sim¬ 
ulations is not far from being realistic. This also agrees with the observations 
of synchrotron and inverse-Compton emission of the Crab nebula, which show 
that the magnetic field is energetically sub-dominant to the population of rel¬ 
ativistic electrons by a factor of 30 dffillas et al.|[l9Mt . 


The above results suggest that the magnetic dissipation inside PWNe is a 
key factor of their dynamics. Combined with the magnetic dissipation in the 
striped zone of the pulsar wind, it allows to reconcile the observations of the 
Crab nebula with the expected high magnetization of such winds, thus finally 
resolving the long-lasting cr-problem. 
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Fig. 18 Depende nce of the total pressure distribution, Igptot, on the imposed symmetry 
( [Forth et al.|2014| . The upper panel shows the pressure distribution in the xz plane of the 
3D simulations and the lower panel in the corresponding 2D run. One sees that the strong 
axial compression observed in 2D simulations is an artifact of the imposed symmetry. 
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2.3 The structure of the highly magnetized region in PWNe 


Even after decay of the alternating fields, the pulsar wind remains highly mag¬ 
netized at high latitudes. The energy and momentum flux in this domain of 
the wind is relatively small therefore the termination shock approaches signif¬ 
icantly closer to the pulsar near the axis than at the equator. This region is 
not quite well resolved in simulations; they just show rapid disruption of the 
flow by the kink instability. On the other hand, the highly magnetized region 
is of special interest because the recently discovered strong, short gamma-ray 
flares from the Crab nebula ( Abdo et al.|2Qllb Tavani et al.|2QH ) are gener¬ 
ally attributed to a rapid magnetic energy release via, e.g., reconnection (see 
Sect. [^, w hich assumes a magnetically dominated region within the nebula. 
Lyubarsky (2012) developed a simple model clarifying the structure of the 


high latitude flow. 

In a highly magnetized flow, only weak shocks could arise therefore the ter¬ 
mination shock at high latitudes is weak; and the postshock flow in this region 
remains radial and relativistic. Lyubarsky (2012) found a simple analytical 
solution for such a relativistic postshock flow. According to this solution, the 
flow initially expands and decelerates but eventually becomes to converge be¬ 
cause the magnetic hoop stress is not counterbalanced either by the poloidal 
field or by the plasma pressure. In the converging flow, magnetic energy is 
converted into the plasma energy therefore the plasma accelerates and heats. 
If the flow remained axisymmetric, it would eventually be focused at the axis, 
the magnetic energy being transferred to the plasma. The focus occurs on the 
axis of the system at the distance from the pulsar ^ where is the open¬ 
ing angle of the highly magnetized part of the wind, a the equatorial radius 
of the termination shock. This point may be identified with the base of the 
observed jet. 


An important point is that in a converging flow, even small perturba¬ 
tions eventually destroy the regular structure. The reason is that if converging 
loops are initially shifted one with respect to another by a displacement much 
less than their radius, the distortion becomes strong when the radius of the 
loops approaches the initial displacement. One has to conclude that when 
the axisymmetric flow is focused into a point at the axis, the magnetic loops 
inevitably come apart close enough to the converging point giving rise to a 
specific turbulence of shrinking magnetic loops. Hence one can expect that 
the energy of the highly magnetized part of the pulsar wind is released in a 
small region close to the converging point; this gives rise to the observed jet. 
Relativistic turbulent motions in highly magnetized plasma imply E ^ B so 
that in the energy release region, particles could be efficiently accelerated ei¬ 
ther via the second order Fermi mechanism or via the magnetic reconnection. 
Therefore the synchrotron gamma-ray emission in the hundreds MeV band, 
both persistent and flaring, could come from a small region at the base of the 
jet. 
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2.4 The unsolved problem: origin of PWN spectra 


One sees that the overall morphology of PWNe is now well understood in the 
scope of MHD models. However, our ignorance of the physical processes giving 
rise to particle acceleration forces us to treat the injection particle spectra in 
PWNe as free parameters ( Del Zanna et al.|2QQ6 Volpi et al.|2QQ8 Olmi et ^ 
2014), and this freedom in interpreting the data limits the level of scrutiny to 


which MHD models can be subjected. The radiation spectrum carries informa¬ 
tion about the particle acceleration processes. Most of the observed radiation 
(from the radio up to a few hundred MeV) is synchrotron emission, with only 
the peak in the very high energy gamma-ray band being attributed to the 
inverse Compton scattering of synchrotron photons off high-energy electrons 
and positrons. The observed strong polarization in the radio, optical and X- 
ray bands is a supporting evidence for the synchrotron origin of the nebula 
emission. 

The synchrotron part of the spectrum may be described as a broken power- 
law. The generic observational feature of PWNe is a flat radio spectrum, Jv oc 
with a between 0 and 0.3, extending in some cases out to the infrared. 
At high frequencies, the spectrum softens, and in the X-ray band, a> 1. Such 
an injection spectrum suggests an unusual acceleration process. The observed 
radio spectrum implies a power-law energy distribution of injected electrons, 
N{E) oc with a shallow slope 1 < k < 1.6. Such an energy distribution 
is remarkable in that most of the particles are found at the low energy end 
of the distribution, whereas particles at the upper end of the distribution 
dominate the energy density of the plasma. Specifically in the Crab Nebula, the 
observed emission spectrum implies that the particles in the energy range from 
^min <100 MeV to ^break ^ 1 T^V are injected into the nebula with a spectral 
slope K = 1.6, so most of the injected energy 5*10^^ erg/s) is carried by TeV 
particles, whereas ^100 times more particles are found at low energies of less 
than 100 MeV. This means that the acceleration process somehow transfers 
most of the total energy of the system to a handful of energetic particles, 
leaving only a small fraction of the energy for the majority of the particles. 
This is not what one would normally expect from the conventional first-order 
Fermi acceleration process, in which the particle flow is randomized at the 
shock and only a fraction of the upstream kinetic energy is deposited in highly 
accelerated particles. 


Lyubarsky (2003b) proposed that the flat energy distribution is formed in 


the course of the particle acceleration by driven reconnection of the alternating 
magnetic field at the pulsar wind termination shock. As a model for such a pro¬ 
cess, Lyubarsky and Liverts (2008) have performed PIC simulations of driven 


magnetic reconnection in a pair plasma. Two stripes of opposite magnetic po¬ 
larity were compressed by means of an external force, which would imitate the 
effect of a shock. They found that driven magnetic reconnection can produce 
flat non-thermal particle spectrum, with k, ^ 1 . Realistic 3D PIC simula¬ 
tions of the shock in a striped wind were performed by [Sironi and Spitkovsky| 
(2011a). They found that the spectrum of accelerated particles depends on the 
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parameter ^ = Xjari^ where A is the stripe wavelength, = mc^rjeB the 
Larmor radius corresponding to the upstream Lorentz factor of the flow, T, 
and the upstream magnetic field, B. It turns out that broad particle spectra 
with flat slopes {1 < k < 2 ) could be formed by the shock-driven reconnection 
only if the above parameter is not less than a few hundreds. In the opposite 
case, the spectrum resembles a Maxwellian distribution. 

Note that aB is in fact the Lorentz factor the particles would achieve if 
the whole spin-down energy were equally distributed between them. Therefore 
the parameter ^ is in fact equal, to within a numerical factor, to the ratio of 
the pulsar light cylinder radius to the Larmor radius acquired by the particles 
when a significant fraction of the Poynting flux is converted to the plasma 
energy. The latter is generally very large because the magnetic field at the 
pulsar wind termination shock is weak therefore ^ could hardly exceed unity. 
The condition for the formation of the non-thermal tail could be achieved only 
if the pulsar wind is overloaded by pairs; then the energy per particle may be 
small enough so that the Larmor radius remains small. At present, there is no 
reason, neither observational not theoretical, to believe in such an necessary 
extraordinary large pair production in pulsars. Therefore the problem of flat 
particle energy spectra in PWNe remains unsolved. 


3 Implications of the Crab flares 


The rapid variability is now a well-est ablished intrinsic property of the Crab 
Nebula in the GeV gamma-ray bancf^ ( Tavani et al.|2Qir Abdo et al.|2Qllb ). 
The flares were not predicted by the models and they generally do not fit 
in the framework of the classical theory of pulsar wind nebulae and particle 
acceleration. This unexpected phenomenon is also a challenge for observers, 
because the Nebula is routinely used as a standard candle for cross-calibrating 
X-ray and gamma-ray instruments. We explain why the flares are so challeng¬ 
ing for the models of the Crab Nebula in Sect. |3.l| and we briefly review some 


of the current attempts to model the flares in Sect. 3.4 (see also the reviews 
by Arons]|2Q12 Biihler and Blandford||2Q14 ). 


3.1 The puzzling features of the flares 


With more than 6 years of data, we know that the mysterious engine at the 
origin of the gamma-ray flares turns on about once or twice a year for about 
a weelj^ Outside of these spectacular events, identified as the “flares”, the 
>100 MeV lightcurve remains apparently restless with continuous small vari¬ 
ations of the flux (Buehler et al. 2012; Striani et al. 2013), as if the engine 


The X -ray flux of the Crab Nebu la is also variable but to a ~ 10% level over a year- 
timescale ( |Wilson-Hodge et al.|[20TT] ). This will not be discussed here, because this phe¬ 
nomenon does not appear to be directly connected to the flares. 

See http://fermi.gsfc.nasa.gov/ssc/data/access/lat/msl_lc/source/Crab_Pulsar 
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Fig. 19 Gamm a-ray light curve abo ve 100 MeV of the April 2011 super-flare, measured by 
the Fermi-LAT ( [Buehler et al.|2012| >. The horizontal dashed blue line shows the “quiescent” 
synchrotron flux > 100 MeV. 


never really switches off. The duration of the flares indicates that the emitting 
region must be surprisingly small compared with the size of the Nebula. For 
a typically week-long episode, the length-scale of the accelerator is of order 
ctflare ^ 10^^ cm, i.c., much smaller than the size of the termination shock 
radius, which is of order 0.1 pc. The brightest events present intra-flare vari¬ 
ability timescales as short as about 8 hrs (see Fig. 19, and Balbo et akpOlT 


Buehler et al.|2Qf2 Mayer et al.|2Q13 ), which put even more severe constraints 


on the size of the particle acceleration site. Consequently, a tiny fraction of the 
Crab Nebula is radiating ^10 times more flux than the entire quiescent Nebula 
in the GeV band. During the April 2011 flare, the gamma-ray flux p eaked at 
about 1% of the spin-down power of the pulsar (Lsd = 5 x 10^^ erg/s) ( [Buehler] 


et al. 2012). This is putting strong constraints on the energetic budget required 


to power the flares. 

The gamma-ray flare spectrum appears at the high-energy end of the qui¬ 
escent synchrotron spectrum, and extends up to about 1 GeV. The flaring 
component is usually attributed to synchrotron radiation emitted by 10^^ eV 
(or PeV) electron-positron pairs in a ^ milli-Gauss magnetic field. Other ra¬ 
diative processes such as inverse Gompton scattering or Bremsstrahlung are 
far too inefficient to cool down the particles over the duration of the fiarj^ 
It was already known that the Grab Nebula accelerates particles up to PeV 
energies (e.g., de Jager and Harding||1992 ). What is new, however, is the evi- 


The fl ux of the inverse-Compton component above > 100 Ge V remains constant during 
the flares ( |H E. S. S. Gollaboration et al.|2Q14[ [Alin et al.|2014| >. 
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dence that such particles are accelerated over such a short timescale (the rise 
and decay timescales of the flares range between 6 hours and few days). In 
fact, the gyration time of the PeV particles is of order the duration of the 
flare themselves. Hence, the particles must be accelerated over a sub-Larmor 
timescale, i.e., the acceleration process must be extremely efficient. Diffuse 
shock-acceleration is not adequate to explain the flares because it operates 
only over multiple gyrations of the particles moving back and forth through the 
shock front. In addi tion, the inferred fla ring particle spectrum can be as hard 
as dN/dj (X 7 “^-^ (Buehler et al. 2012) which is inconsistent with the steep 
power-law (i.e., of index ^ — 2 ) usually expected in diffuse shock-acceleration 
(e.g., [Blandford and Eichler||1987| Sironi and Spitkovsky||2QQ9a ). 


3.2 Multiwavelength follow-up of Crab PWN flares 


During the gamma-ray flares the Crab was observed in radio, infrared, optical. 
X-ray and TeV energies, but no substantial variation in the flux emission at 
these wavelengths was measured. After the discovery of the gamma-ray flares, 
the Chandra X-ray observatory started to observe the Crab approximately 
every month. Five observations were carried out during the major gamma-ray 
flare of April 2011. The bright Anvil region (see Figure and several other 
regions (that are known to be active) exhibit time variability during the flaring 
activity time. Nevertheless, despite these hints in the X-ray data, there is no 
evidence for statistical significant variations associated with the flare |Weisskopf| 


et al. (2013). The near-IR observations performed by Keck’s NIRC 2 revealed 


that the inner knot (knot- 1 ) was slightly brighter when compared to previous 
observations. Indeed knot-1, which is the brightest feature from the Nebula 
in the near infrared energy band, was reported to show flux variation at this 
wavelength of the order of 35%. However this variation is well within the 
range typically observed from this region. The radio observations performed 
with VLA did not reveal anything interesting. No other point source, a part 
from the pulsar, was found. Therefore, no “smoking gun” has been identified 
from the X-ray, near-IR and radio observations. 

In terms of X-ray and optical counterparts of the Crab flares, we can de¬ 
scribe two scenarios: 


1 . A simultaneous brightening of X-rays and optical locations associated with 
the gamma-ray flares. This scenario would favor “shock-driven” power-law 
models of particle acceleration; 

2 . A delayed response of the optical/X-ray emission, with a timescale that 
depends on the radiative cooling properties of the accelerated particle pop¬ 
ulation. This scenario would favor an extremely efficient acceleration mech¬ 
anism, likely to saturate the particle energy to a maximum value, and that 
could be modeled with a quasi-monoenergetic particle distribution; 

The absence of a strongly enhanced X-ray and/or optical location in the 
Nebula in coincidence with the gamma-ray flares tends to exclude the first 















40 


Oleg Kargaltsev et al. 


scenario. However in Weisskopf et al. (2013) the possibility to model the photon 
spectrum of the flare of April 2011 with a power-law with index F ^ 1.3 
connecting gamma rays and X rays was investigated and could not be ruled 
out. 

Because of the poor angular resolution of gamma-ray telescopes, localizing 
the emission site of the Crab flares is a big challenge. Several regions can be 
considered as candidate for the acceleration and emitting region. Among them, 
we can identify 3 particularly interesting possibilities: 

1. instabilities in the Anvil at the South-East jet base. Variability in this 
region was detected by optical and X-ray observations both during the 
September 2010 and April 2011 gamma-ray flares ( Tavani et al.|20lT Weis¬ 


skopf et al. 


2013). MHD simulations of the Crab South-East jet (Mignone 


et al. 2013) revealed substantial jet deviation and magnetic dissipation. 


The jet, fed by highly magnetized and relativistic plasma (a ^ 1 — 10, 


Mignone et al. 2013), could be a region of magnetic reconnection. Kink 


instability in the jet could trigger magnetic reconnection and consequent 
particle acceleration. 

2. tearing mode and reconnection on the termination shock “ring”. This re¬ 
gion is known to be highly variable. The three most variable spots during 
the April 2011 flare are located along the ring. The highly variable wisps 
are observed to originate from this region. Recent simulations ( [Forth et al.| 


2014) recently proposed that efficient magnetic reconnection could take 


place in the Nebula right after the termination shock; 
variation in the observed emission from knot-1. Variability up to ^ 20—30% 
is known. Its variability is interpreted as a variable Doppler factor but 
no indication for optical variability has been observed so far. Simulations 


(Camus et al. 2009 Komissarov et al. 2007) show that at 100 MeV the 


inner knot is the brightest feature of the Nebula, and that the magnetic 
field can be up to 10 times larger than the average Nebular magnetic field. 


3.3 Proposed models of the flares 


It is usually expected that synchrotron spectrum observed in astrophysical 


sources cuts off far below 100 MeV (Blumenthal and Gould 1970 Guilbert 


et al. 1983). This limit is given by the balance of two antagonist forces acting on 


the particles: (i) the accelerating electric force, and (ii) the radiation-reaction- 
force opposite to the particle’s motion due to the emission of synchrotron 
photons. Hence, there is a maximum energy limited by synchrotron losses 


rather than by the size of the accelerator (e.g., Aharonian et al. 2002 Medvedev 


2003). Then, one can show that the maximum (critical) synchrotron photon 


energy should be 


sync _ 
'-max 


9meC^ 

4(ap 



160 X 



MeV, 
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Gamma-ray flare spectra 



Fig. 20 Measured gamma-ray spectra of the Crab Nebula in the 1 MeV-10 GeV band. 
Black data-points show the “quiescent” spectra of the Nebula. The green, magenta, and 
blue points are respectively the > 100 MeV spectra duri ng the April 2011, September 2010 
and February 2009 flares measured by the Fermi-LAT ( [Blihler and Blandford|[2014| ). The 
red solid line is the spectra obtained with 3D PIC simulations with radiation reaction of 
a reconnecting current sheet in the Crab Nebula ( [Cerutti et al.|2014| ). The model assumes 
that the Nebula is 2 kpc away. 


where tUq is the rest mass of the electron, ~ 1/137 is the fine structure 
constant, E the electric field, and B± the magnetic field perpendicular to the 
particle’s direction of motion. A particle accelerated above the radiation reac¬ 
tion limit would radiate away most of its energy within a sub-Larmor cycle. In 
most cases, ideal MHD applies, E so the synchrotron spectrum should 

cut ^100 MeV. This rule of thumb generally applies well to astrophysical 
sources, and in particular to the quiescent Crab Nebula where the synchrotron 
spectrum turns over at about the 100 MeV limit ( de Jager et al.p996 Abdo 


et al.||2010a ). However, the flaring emission is systematically extending signif¬ 


icantly above 160 MeV, up to about 1 GeV (see Fig. [^. Unless the emitting 
region moves at highly relativistic speeds 0.9 c), it implies that E ^ 
suggesting that a non-ideal MHD process may be at work. This is, once again, 
difficult to explain with classical models of particle acceleration. 

All of the estimates derived above from observations are quite conserva¬ 
tive, in a sense that we ignored the effect of beaming (geometrical or rela¬ 
tivistic), spatial and/or temporal inhomogeneities. Every model proposed so 
far is taking advantage of one or more of these effects to alleviate the tight 
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constraints imposed by the flares. For instance, one feature commonly invoked 
in models is a strong inhomogeneity of the flaring region, in particular in the 
magnetic field structure. Bykov et al. (2012) proposed that the flares occur 


around the equatorial belt where the annihilation of the striped pulsar wind 


at the shock accelerates particles (Lyubarsky 2003b Petri and Lyubarsky 


2007 Sironi and Spitkovsky 2011a) and, crucially for this model, generates 


magnetic turbulence. They find that a concentration of fluctuating magnetic 
field can generate an intermittent, strongly polarized gamma-ray signal that is 
most pronounced at the high-energy end of the synchrotron spectrum. In this 
model, synchrotron photons above 100 MeV can be emitted if the magnetic 
field varies over a timescale shorter than the synchrotron cooling time of the 
particles which is determined only on the RMS value of the field. The observed 
gamma-ray variability would then be given by the statistical properties of the 
magnetic fluctuations (lifetime, amplitude). 

The model by Teraki and Takahara| ( |2013 ) also relies on a highly inho¬ 
mogeneous turbulent flow, but in which the coherence length of the magnetic 
field, Ab, is extremely short compared with the cooling length of the particles 
and even compared with the formation length of the synchrotron photons, 
Async = IeB > Ab, □ In this regime, the particles emit the so-called 
“jitter” radiation rather than the classical synchrotron radiation (see, e.g., 
Medvedev||2QQQ ). The cooling rate of the particles is identical to synchrotron 
but the emitted spectra can differ significantly. In particular, the critical pho¬ 
ton energy increases by a factor Cjitt/csync ^ Agync/AB ^ 1 in the jitter regime, 
and hence > 100 MeV gamma rays could be emitted even by particles below 
the radiation reaction limit. In addition, jitter radiation can produce a harder 
spectrum than synchrotron, specifically, oi v instead of oc for a 
mono-energetic population of particles. However, we note that jitter radiation 
is not needed to explain the flare spectra. So far, observations are fully consis¬ 


tent with synchrotron radiation (Buehler et al. 2012 Weisskopf et al. 2013) 


but this could be tested in the future. 

A more natural way to explain the Crab flares is to invoke a relativistic 
bulk motion of the flaring region with a modest Lorentz factor F ^ 2. In¬ 
deed, the relativistic motion of the source can boost ^ 100 MeV synchrotron 
photons emitted in the co-moving frame by a ^ T factor above the radiation 


reaction limit (Yuan et al. 

2011 Bednarek and Idee 

2011 

Komissarov and 

Lyutikov 

2011 

Lyutikov et al. 

2012 Clausen-Brown and Lyutikov 2012). A 


Doppler boost would also relax the tight constraints on the size of the region 
and the duration of the flare in the co-moving frame, and beam the emission 
in the frame of the observer which would also help at reducing the energetic 
constraint^^ Lastly, the Doppler beaming would explain the observed correla¬ 
tion between the gamma-ray flux and the cut-off energy in the spectra during 
the April 2011 flare ( Buehler et al.|2Ql2 ). Although this simple solution solves 
many problems at once, there is still no definite evidence that such a relativis- 


21 

22 


This might be challenging to achieve in relativistic shocks (Sironi and Spitkovsky 


2009b). 


This solution is often proposed to account for the ultra-rapid gamma-ray flares in blazars. 
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tic flow exists in the Crab Nebula. Observations show only mildly relativistic 
flows with proper motion of order half the speed of light (e.g., Hester [TqQ 8[ ) . 
However, in principle, highly relativistic flows could emerge in the polar regions 
of the Nebula because the relativistic shock is oblique and magnetized (weak 
shock) at high latitudes (e.g., Lyubarsky 2012 Komissarov 2013). In partic¬ 
ular, Komissarov and Lyutikov (2011) argued that the flares may originate 


from the emission from the oblique shock Doppler-boosted towards the ob¬ 
server, that they associate with the well-known bright compact structure near 
the pulsar (the so-called “inner knot”). Unfortunately, the brightness of the 
knot does not show any variations correlated with the gamma-ray flares (see 
ITavani et al.||201l| [Lobanov et al.||2QTT| jWeisskopf et al.|[2Q13| and references 
therein) contrary to what the model predicts. 


Alternatively, Uzdensky et al. (2011) proposed that magnetic reconnection 


could accelerate particles well above the radiation reaction limit, and hence 
result in the emission of > 100 MeV synchrotron radiation. Indeed, as pointed 
out by Kirk (2004), inside a reconnection layer the magnetic field is small and 


even vanishes at its center while the reconnection electric field is maximum, 
i.e., we are in the situation where E ^ B. Thus, in principle, a particle trapped 
deep inside the reconnection layer could be linearly accelerated by the electric 
field to arbitrary high energies with little synchrotron losses. The maximum 
energy of the particle would be limited only by the length of the reconnection 
layer L, i.e., fmax ^ eEL = e^recBoL where Pj.ec = EIBq is the dimensionless 
reconnection rate, E is the reconnection electric field and Bq is the recon¬ 
necting magnetic field. Using a test-particle approach with prescribed static 


fields, Contopoulos (2007) and Cerutti et al. (2012a) showed that the high- 


energy particles are naturally trapped and confined deeply inside the layer 
where they follow the relativistic analog of Speiser orbits (Speiser 1965 Uz- 
densky et al. ]MT] ). This scenario was successfully tested using state-of-the-art 


particle-in-cell (PIC) simulations of 2D and 3D relativistic reconnection with 
guide field and, more importantly, with the radiation reaction force (Cerutti 
et al.|2013 , 2014). Furthermore, these studies revealed that a natural outcome 


of relativistic reconnection is the strong beaming and bunching of the ener¬ 
getic particles (see also Cerutti et al.||2012b| ). The combination of both effects 


results in several bright ultra-rapid synchrotron flares above 100 MeV consis¬ 
tent with the observed intra-flare variability (^ 6 hrs) if, by chance, the beam 
crosses the line of sight of a distant observer. Thanks to beaming, this model 
can also explain the overall energetics of the flares (Fig. [^, the flux/cut-off 
energy correlation and, at least qualitatively, the hard particle spectrunj^ 
The reconnection scenario works best in a highly magnetized flow (i.e., 
cr ^ 1) which may be hard to find in the nebula, except in the polar regions 


and in the jets. As predicted by I 

Begelman 

(1998 

) and as recently shown by 

Mizuno et al. (2011); Forth et al 

. (2014); 

Mignone et al. (2013 

), the polar 


regions and the jets are subject to kink instabilities which results in important 


Recent studies sh ow that reconnection produces hard particle spectra dN / d'y oc 
7 “^, 7 “^-^ for (7 ^ 1 ( Sironi and Spitkovsky|2oT4 Guo et al.|2014 Werner et al.|2Q14 >. 
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magnetic dissipation (see Sect. 2.2), and may ultimately power the Crab flares 


( Cerutti et al.|2Ql^ Lyubarsky|2Q12 Komissarov|2Q13 Mignone et al.|2QT3 ) 


3.4 Comparison with GeV flares in the PSR B1259-63/LS 2883 binary 


The gamma-ray binary PSR B1259—63 that contains a 48-ms pulsar in a 
3.4-yrs eccentric orbit around an O star (see Sect. 1.4), is a well-known TeV 
gamma-ray emitter (Aharonian et al. ]M^|H.E.S.S. Collaboration et al.|201^ . 
The very-high energy radiation is often interpreted as inverse Compton scat¬ 
tering of the UV stellar photons on relativistic pairs accelerated near the shock 


front between the pulsar wind and the stellar wind (e.g., Tavani and Arons 


1997 Kirk et al. 1999; Khangulyan et al. 2007). The model predicts a maximum 


of GeV gamma-ray emission close to periastron which was indeed observed for 


the first time by the Fermi-LAT during the 2010 periastron passage (Tam et al 
2011 Abdo et al.|2011a ). A few weeks after the peak of the periastron emission 
faded away, and against all expectations, a bright flare appeared in the Fermi 
data. The flare is about 10-times brighter than the predicted emission at pe¬ 
riastron, which represents a gamma-ray luminosity comparable to the pulsar 
spin-down power. A radiative efficiency close to 100% is in principle achievable 
with inverse Gompton scattering, but the density of stellar photons is far too 
low to explain the flux at these phases, unless there are extra sources of radia¬ 
tion close to the pulsar (see, e.g., Khangulyan et ar]|2012 Dubus and Cerutti 


2013), or significant Doppler boosting of the emission towards the observer 
( Dubus et al.||2010 Kong et al.|2012 ). The flare peaks at 300 MeV and is seen 
only in the GeV band, which suggests that the particle energy distribution 
must be very narrow. These properties remind us of the Crab-flare events (see 
Sect. 3.1). Although both flares share similar properties, they have also im¬ 


portant differences. The week-long flares in the Crab reach at most 1% of the 
Crab pulsar spin-down power, whereas the highest day-aver age flux reaches 
nearly 100% of the pulsar spin-down flux during flaring period associated with 
the second disk passage. Unlike Crab PWN the flares in B1259-63 appear to 
be periodic, i.e. they occurred at similar binary phases (close to the second 
disk passage) during the past two binary cycles (Tam et al.||2015 Caliandro 


While it may be that the underlying nature of flares in both cases is the 
synchrotron radiation associated with reconnection in the magnetized rela¬ 
tivistic plasma, the details should differ. In the Crab PWN, the reconnection 
can be driven by the growth of instabilities or other random process, while in 
B1259-63, it can be driven by the magnetic field distortion and compression 
caused by the pulsar passage through the excretion disk. It is possible that the 
reconnection happens in the tail of the PWN after the pulsar passage through 
the disk which could explain the delay between the GeV flare and the peak of 


X-ray flux (Tam et al. 

2015 Galiandro et al.||2015 

). An alternative scenario. 

considered by Khangulyan et al. 

(2012 

), where the flare is due to the IG radia- 


tion associated with the increase in volume occupied by the unshocked pulsar 
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wind when the excretion disk is strongly perturbed, also remains a possibility. 
One may be able to differentiate between the two scenarios once the 7 -ray 
variability timescales and the state of the disk during the pulsar passage are 
better probed by the observations. 
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